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AFTOX  4.0  -  The  Air  Force  Toxic  Chemical 
Dispersion  Model  -  A  User's  Guido 


1.  INTRODUCTION 


The  AFTOX  model  Is  an  interactive  Gaussian  puff/plume  dispersion  model  designed  to  run 
on  a  desk  top  computer.  The  program,  written  in  GW-BASIC  language.  Is  user-oriented  and 
written  in  a  conversational  mode  to  allow  its  operation  with  only  a  minimum  knowledge  of 
the  computer  or  the  model.  The  model  is  designed  to  handle  continuous  and  Instantaneous, 
liquid  and  gas.  surface  and  elevated  releases  from  a  point  or  area  source.  The  model  will  1) 
plot  concentration  contours  and  toxic  corridor,  2)  compute  concentration  at  a  specified  point 
and  time,  and  3)  compute  maximum  concentration  at  a  given  elevation  and  time.  It  contains  a 
library  of  129  chemicals  but  may  be  run  for  other  chemicals  as  well.  The  model  also  predicts 
the  dispersion  of  a  continuous  buoyant  plume  from  a  stack. 

The  model  has  many  unique  features,  including  the  computation  of  a  continuous  stability 
parameter,  the  inclusion  of  the  concentration  averaging  time  and  surface  roughness,  the 
ability  to  print  now  or  save/prtnt  later,  the  ability  to  correct  input  data  without  restarting, 
easy  access  to  the  chemical  data  file  for  adding,  deleting  and  changing  data,  and  a  plot  of  the 
plume  and  90  percent  hazard  area  with  automatic  scaling. 

The  AFTOX  model  was  developed  as  a  result  of  a  need  by  the  U.S.  Air  Force  to  update  its 
toxic  corridor  prediction  capability.  The  Air  Force  had  been  using  the  Ocean  Breeze/Dry  Gulch 
model,  which  was  an  empirical  model  derived  In  the  mid-sixties  from  two  9ets  of  diffusion 
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experiments:  Ocean  Breeze,  conducted  at  Cape  Canaveral,  Florida  and  Dry  Gulch,  conducted  at 
Vandenberg  AFB,  California. 

The  requirements  for  the  new  model  were:  1)  that  it  be  able  to  run  on  a  microcomputer,  2) 
that  it  be  very  user-friendly.  3)  that  the  Input  data  be  limited  to  data  that  would  normally  be 
available  to  the  operator,  and  4)  that  the  computation  time  be  kept  to  a  minimum. 

To  meet  these  criteria,  it  was  decided  that  a  Gaussian  puff  dispersion  model  would  be  most 
appropriate.  As  a  result,  AFTOX  was  developed  and  evaluated  against  243  diffusion  tests  from 
the  Prairie  Grass, 1,2,3  Ocean  Breeze  and  Dry  Gulch,4,6  and  Green  Glow6  experiments.  In  1988, 
the  model  was  officially  endorsed  by  the  Air  Weather  Service,  and  Is  now  used  extensively 
throughout  the  Air  Force. 

AFTOX  4.0  is  an  updated  version  of  the  original  model.7  Many  of  the  changes  made  to  the 
original  version  were  a  result  of  suggestions  made  by  the  users.  Many  of  these  changes  are 
designed  to  reduce  the  chances  of  Input  error  and  to  provide  more  default  values.  The  most 
significant  technical  changes  are  1)  an  expansion  of  the  chemical  data  base  from  76  to  129 
chemicals:  2)  the  inclusion  of  the  more  universal  Antoine  equation  for  computing  the 
saturated  vapor  pressure,  thus  making  it  possible  to  include  more  chemicals  In  the  data  base; 

3)  a  new  evaporation  model  for  the  hypergols  and  nitrogen  tetroxide:  4)  a  slight  adjustment  In 
the  vertical  dispersion  coefficients  to  bring  them  more  In  line  with  the  Pasqulll-Glfford8 
dispersion  coefficients:  and  5)  the  Inclusion  of  the  90  percent  hazard  area  for  all  spills  and  not 
Just  for  continuous  releases  as  in  the  original  version. 

Section  2  of  this  report  presents  a  technical  description  of  the  model.  This  section  Is 
similar  to  the  same  section  In  the  previous  user’s  guide.7  but  with  some  Important  differences. 
The  entire  technical  discussion  is  repeated  here  for  completeness.  Section  3  Is  a  user's  guide  to 
the  operation  of  the  model. 


1  Barad,  M.L..  Ed.  (1958)  Project  Prairie  Crass,  A  Field  Program  in  Diffusion.  Vol.  I.  AFCRC-TR- 
58-235(1).  AD  152572. 

2  Barad,  M  L..  Ed.  (1958)  Project  Pratrle  Grass,  A  Field  Program  in  DiJJusioix.  Vol.  II,  AFCRC- 
TR-58-235(II),  AD  152573. 

3  Haugen.  D.A.,  Ed.  (1959)  Project  Prairie  Grass,  A  Field  Program  in  Diffusion,  Vol.  III. 

AFCRC -TR-58-235(III),  AD  217076. 

4  Haugen.  D.A  and  Fuquay,  J.J..  Eds.  (1963)  The  Ocean  Breeze  arid  Dry  Gulch  Diffusion 
Programs,  Vol.  1,  AFCRL-63-791U).  AD  428436. 

5  Haugen.  D.A.  and  Taylor,  J.H..  Ed9.  (1963)  The  Ocean  Breeze  and  Dry  Gulch  Diffusion 
Programs.  Vol.  2,  AFCRL-63-79 1(11).  AD  427687. 

6  Barad.  M  L.  and  Fuquay,  J.J.,  Eds.  (1962)  The  Green  Glow  Diffusion  Program,  Vol.  1  &  II, 
AFCRL-62-25(I  &  III. 

7  Kunkel,  B_A.  (1988)  User's  Guide  for  the  Air  Force  Toxic  Chemical  Dispersion  Model 
(AFTOX),  AFGL-TR-8 8-0009.  ADA  199096. 

8  Turner,  D.B.  (19/0)  Workbook  of  Atmospheric  Dispersion  Estimates.  Public  Health  Service 
Publication  No.  999-AP-26  (PB  191  482),  U.S.  Dept,  of  Health,  Education  and  Welfare. 
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2.  THEORY 


2.1  Basic  Diffusion  Equations 

The  Gaussian  puff  model  uses  an  equation  to  describe  the  dispersion  of  a  puff  with  time. 
The  equation  assumes  that  the  material  is  conserved  during  transport  and  diffusion,  that  Is. 
there  is  no  decay  or  deposition.  It  further  assumes  that  the  distribution  of  concentration 
within  the  puff  is  Gaussian. 

The  Gaussian  puff  equation  can  be  written 


G  (x.  y.  z.  t  -  t')  = 


9(0 


o  3/2 
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where  G  is  the  concentration  in  the  puff  at  a  given  point  (x,y,z)  and  time  (t  -  t').  Q  is  the  total 
mass  In  the  puff.  The  diffusion  parameters.  ox,  oy,  oz  are  the  standard  deviation  of  the 
material  concentration  in  the  x.  y.  and  z  directions.  It  is  assumed  that  ax  =*  oy.  thus  producing 
a  circular  horizontal  puff  cross- section.  The  variable  t  represents  the  total  elapsed  time  since 
the  spill,  and  the  t'  is  the  time  of  emission  of  the  puff.  Thus,  (t  -  t')  is  the  travel  time  or 
elapsed  time  since  the  puff  emission;  u  is  the  wind  speed  at  10  m;  and  H  Is  the  height  of  the 
source. 

When  there  is  an  Inversion,  Eq.  (1)  is  modified  by  adding  the  following  expression  to  the 
last  two  terms  In  Eq.  (1) 
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where  L  Is  the  mixing  layer  height  and  N  is  the  number  of  iterations.  The  computer  program 
considers  the  series  converged  when  the  above  expression  has  a  value  less  than  0.01  for  a  given 
N.  If  there  is  no  Inversion,  the  expression  is  Ignored  in  the  program. 
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When  the  series  requires  a  large  number  of  Iterations  to  converge  (>200)  the  program 
assumes  a  uniformly  distributed  nUme  between  the  earth's  surface  and  the  Inversion  height.  If 
this  Is  the  case,  then  the  concentration  at  a  given  point  and  time  In  the  puff  Is  described  by  the 
expression: 


G  (x,  y.  z.  t  -  t')  = 


0(f) 


2no  o  L 
x  y 


exp 
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The  concentration  at  a  point  In  space  at  a  given  time  depends  on  the  number  of  nearby  puffs, 
their  size,  and  the  amount  of  material  In  each  puff.  The  sum  effect  of  all  these  puffs  Is  given 
by  summing  over  all  emission  times: 


oo 


G  (x.  y.  z.  t)  =  £  G  (*•  y  z  t~  f  )■ 
t*  o 


(4) 


For  an  instantaneous  gas  release  there  Is  only  one  emission  time  and  one  puff.  Therefore, 
a  summation  Is  not  necessary.  However,  for  a  continuous  spill,  or  a  spill  of  flnile  duration, 
the  summation  Is  performed  over  puffs  whose  centers  are  located  within  four  standard 
deviations  of  the  puff  concentrations  upwind  and  downwind  from  the  location  of  interest.  It  is 
assumed  that  concentrations  from  puffs  further  than  four  standard  deviations  contribute  little 
to  the  concentration  at  the  specified  location. 

In  a  spll)  of  finite  duration,  the  model  assumes  20  puffs/min.  out  to  300  m  from  the  source, 
4  puffs/ min.  from  300  m  to  3  km  and  then  3  puffs/mln.  beyond  3  km  for  winds  less  than  or 
equal  to  4  m/sec.  For  winds  greater  than  4  m/sec,  the  number  of  puffs  per  minute  (N)  Is  equal 
to 


N  =  INT  (N,  u/4). 


(5) 


where  N,  Is  20,  4.  or  3  depending  on  the  distance  from  the  source.  The  reason  for  Increasing 
the  frequency  of  puffs  Is  that  higher  winds  cause  greater  distances  between  pufls.  thus  reducing 
the  number  of  overlapping  puffs.  By  Increasing  the  frequency  of  puffs  according  to  Eq,  (5), 
approximately  the  9ame  number  of  overlapping  puffs  Is  retained  no  matter  how  strong  the 
wind. 

Under  steady-state,  non-inversion  conditions,  keeping  track  of  the  Individual  puff9  is  not 
necessary.  A  simple  Gau99lan  plume  model  will  suffice.  In  this  situation,  the  model  uses  the 
following  Gaussian  plume  equation  in  place  of  Eq.  (1): 
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2.2  Atmospheric  Stability 

The  atmospheric  stability  and  the  distance  from  the  source  affect  the  dispersion 
parameters  (oy,  o2).  In  diffusion  modeling,  the  atmospheric  stability  is  often  defined  by  the 

Pasqulil  stability  categories,  which  range  from  category  A  for  a  very  unstable  atmosphere  to  F 
for  a  very  stable  atmosphere.  In  AFTOX,  a  continuous  stability  parameter  ranging  from  0.5  to 
6  is  used  in  place  of  the  discrete  stability  categories.  This  prevents  sharp  changes  in  the 
hazard  distance  when  going  from  one  stability  category  to  another;  this  can  happen  with  a 
slight  change  in  wind  speed,  solar  angle,  or  cloud  cover.  The  relationship  between  Pasquill 
stability  category  and  the  continuous  stability  parameter  is  shown  in  Table  1. 


Table  1.  Relation  Between  the  Pasquill  Stability  Categories.  SC.  and  the  Continuous 
Stability  Parameter.  SP 

jmb— bw awapacBB—ne 


SC  A 

B 

C 

D 

E 

F 

SP  0-1 

1-2 

2-3 

3-4 

4-5 

5-6 

AFTOX  employs  one  of  two  methods  to  define  the  stability  parameter:  (1)  uses  the  wind 
speed  and  solar  insolation,  or  (2)  uses  the  standard  deviation  of  the  wind  direction  to  define 
the  stability  parameter.  In  the  first  method,  the  stability  is  determined  using  Golder's 
nomogram9  where  the  Monln-Obukhov  length  (L)  and  surface  roughness  are  related  to  the 
Pasquill  stability  categories.  Since  L  is  not  measured,  it  must  be  calculated,  and  because  L  is  a 
function  of  the  friction  velocity  (u.)  and  sensible  heat  flux  (H),  these  latter  two  parameters 
must  also  be  calculated.  In  the  second  method,  the  stability  is  determined  using  the  Modified 
Sigma  Theta  (MST)  approach  Introduced  by  Mitchell  and  Timbre10  and  presented  in  more 
detail  in  a  later  report  by  Mitchell.1 1 


9  Golder.  D.  (1972)  Relations  between  stability  parameters  in  the  surface  layer,  Boundary 
Layer  Met.  3:46-58. 

10  Mitchell.  A.E..  Jr.  and  Timbre.  K  (1979)  Atmospheric  stability  class  from  horizontal  wind 
fluctuation,  paper  79-29.2,  Air  Pollution  Control  Association  Annual  Meeting,  Cincinnati.  OH. 

11  Mitchell.  A.E.,  Jr.  (1982)  A  comparison  of  short-term  dispersion  estimates  resulting  from 
various  atmospheric  stability  classification  methods,  Atmos.  Environ.  16:765-773. 
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2.2.1  CALCULATING  u.  AND  L 


To  determine  u.  and  hence  L,  an  Iterative  procedure  similar  to  that  described  by  Koo  et  al12 
Is  applied  to  the  surface  layer  wind  profile  equations  below,  as  summarized  by  Ragland  and 
Dennis. 13  for  neutral,  unstable,  and  stable  conditions. 


Neutral:  fen.  =  In  -2-  ; 

^  Zo 


(7a) 


Unstable: 


2  (tan  1  x-  tan-1  x  )  +  In  — — — 
v  oJ  x  +  1 


In 

x  +  1 

O 


(7b) 


where  x  =  (1-15Z/L)1/4  and  Xq=  (1-15Z0/L)1/4,  and 


Stable: 


kli  =  In  -Z-  +  5.2  a, 

ZC 


(7c) 


where  u  is  the  measured  wind  speed,  u,  is  the  friction  velocity,  k  is  the  von  Karman  constant,  Z 
Is  the  anemometer  height,  ZQ  Is  the  roughness  length  at  the  wind  measuremenl  site,  L  Is  the 
Monin-Obukhov  length,  and  a  =  Z/L  when  Z  <  L  and  a  =  1  when  Z  >  L. 

The  Monin-Obukhov  length  L  Is  defined  as 


PC  Tu.3 
L  =  £ 

kgH 


(8) 


where  p  Is  the  ambient  air  density,  Cp  Is  the  specific  heat  of  air,  T  is  the  air  temperature,  k  Is 
the  von  Karman  constant,  g  Is  the  gravitational  constant,  and  H  is  the  sensible  heat  flux. 
Equation  (8)  may  be  rewritten 


12  Koo,  E.,  Chang,  W.L.,  and  Lee,  B.Y.  (1984)  An  alternative  stability  classification  scheme 
for  air  quality  modeling,  Fourth  Joint  Conference  on  Applications  of  Air  Pollution 
Meteorology,  Portland,  OR,  95-98. 

13  Ragland.  K.W.  and  Dennis,  B.L.  (1975)  Point  source  atmospheric  diffusion  model  with 
variable  wind  and  dlffusivlty  profiles,  Atmos.  Environ.,  9:175-189. 
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(9) 


L  *  - 


Pu? 


0.01 12H 


where  L  Is  lr  m,  p  Is  the  atmospheric  pressure  In  mb,  u.  is  in  m  sec'1  and  H  Is  in  w  m'2.  The 

method  used  to  compute  the  heat  flux  Is  described  in  Section  2.2.2. 

To  start  the  Iteration,  an  estimate  of  u.  is  required.  That  estimate  is  calculated  from 

Eq.  (10): 


_1 

12 


0.8H  \ 
H  +  1000/ 


(10) 


where  the  wind  speed  at  10  m  (u)  is  estimated  from  the  measured  wind  (uz)  using  the  power  law 
with  an  exponent  of  0.2: 


u  o  u 

t 


(11) 


The  Iteration  continues  until  successive  values  of  u./u  differ  by  less  than  0.001.  L  is  then 
calculated  from  this  value  of  u..  Once  u.  is  derived,  the  wind  speed  at  10  m  Is  determined  and 
used  in  the  Gaussian  puff  equation  lEq.  (1)|  or  plume  equation  (Eq.(6)]. 


2.2.2  CALCULATING  SENSIBLE  HEAT  FLUX 

The  sensible  heat  flux  (H)  Is  needed  for  calculating  the  Monln-Obukhov  length  and  for 
detei  mining  the  appropriate  wind  profile  equation.  A  positive  heat  flux  greater  than  1  W  m'2 
Is  unstable  lEq.  (7b)),  while  a  negative  heat  flux  less  them  - 1  W  m'2  Is  stable  (Eq.  (7c)).  A  heat 
flux  between  1  and  -1  W  m  2  is  neutral  (Eq.  (7a)J.  The  heat  flux  Is  also  used  to  determine  an 
initial  estimate  of  u».  The  sensible  heat  flux  Is  calculated  for  daytime  conditions  using  a 

method  described  by  Holtslag  and  Van  Ulden.14 

The  incoming  solar  radiation  at  ground  level  (SR)  is  very  much  a  function  of  the  solar 
elevation  angle  (<fr).  Under  clear  skies.  It  can  be  determined  from  the  empirical  formula: 


SRC  =  at  sin  +  ♦  a^. 


(12) 


14  Holtslag,  A. A  M.  and  Van  Ulden,  A.P,  (1983)  A  simple  scheme  for  daytime  estimates  of 
the  surface  fluxes  from  routine  weather  data.  J.  Climate  AppL  MeteoroL,  22:517-529. 
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The  turbidity  coefficients,  at  and  a^  describe  the  average  atmospheric  attenuation  of  K0  by 
water  vapor  and  dust  for  a  given  site.  These  coefficients  will  vary  from  one  location  to 
another  and  from  one  time  period  to  another  because  of  the  variation  In  the  turbidity  of  the 
atmosphere.  In  this  model,  aj  =  990  W  m'2  and  a^  =  -30  W  m'2  were  used  since  they  represent  a 

reasonable  average  of  coefficients  determined  at  various  locations.  The  method  of 
determining  the  solar  elevation  angle  ($)  Is  described  in  Section  2.2.3. 

The  presence  of  clouds  reduces  the  Incoming  solar  radiation.  Kasten  &  Czeplak15  proposed 
the  following: 


SR  =  SR0(1  +b1Nb2), 


(13) 


where  N  is  the  fraction  of  the  sky  covered  with  clouds  and  bj  and  b2  are  empirical  coefficients. 
From  ten  years  of  observation  at  Hamburg.  Kasten  and  Czeplak  obtained  b,  =  -0.75  and  b2  =  3.4 
on  the  average.  Using  these  coefficients  Implies  that  the  solar  radiation  reaching  the  ground 
under  overcast  conditions  Is  25  percent  of  that  reaching  the  ground  under  clear  skies.  Kasten 
and  Czeplak  also  show  that  the  transmlttances  (TR)  of  cirrus,  altus.  cumulus,  stratus,  and 
ffimbostratus  are  0.61.  0.27.  0.25.  0.18.  and  0. 16.  respectively.  Therefore,  to  take  Into  account 
the  type  of  cloud.  F.q.  (13)  has  been  modified  to  read: 

SR  =  SRq  (1  -  (1  -  TR)  Nb2).  (14) 

In  the  model,  the  cloud  types  have  been  condensed  Into  three  groups  representing  high, 
medium,  and  low  clouds.  The  cloud  groups  and  their  transmlttances  are  shown  in  Table  2. 


Table  2.  Transmlttances  ior  High,  Medium  and  Low  Clouds 


HIGH 

Cl.  Cc,  Cs 

0.61 

MEDIUM 

Ac.  As,  Sc,  Cu 

0.26 

LOW 

St,  Ns,  Fog 

0.17 

To  determine  the  net  radiation,  the  following  equation  from  Holtslag  and  Van  Ulden  Is 

used: 


15  Kasten,  F.  and  Czeplak,  G.  (1980)  Solar  and  terrestrial  radiation  dependent  on  the 
amount  and  type  of  cloud.  Solar  Energy,  24:177-189. 
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(15) 


Q 


(1  -iOSR  +  CjT^-cT4  +  C2N 

1  +  c„ 


where  r  Is  the  albedo  of  the  surface,  equal  to  0.25,  except  when  there  Is  a  snow  cover,  when  it  Is 
set  at  0.75.  T  Is  the  air  temperature  In  Kelvin,  a  Is  the  Stefan-Boltzmann  constant,  equal  to 
5.67  x  10'8  W  m'2  K‘4.  and  the  empirical  constants  Ct  and  C2  are  equal  to  5.31  x  10  13  W  m'2  K  0 
and  60  W  m'2.  respectively.  The  surface  heating  coefficient  C3  varies  from  0.12  for  a  wet 
surface  to  0.38  for  dry,  bare  soil.  In  the  model.  C3  Is  equal  to  0.12  for  a  wet  surface  and  0.25 
for  a  dry  surface. 

The  surface  energy  budget  relates  the  net  radiation  Q  to  the  various  heat  fluxes  at  the 
earth's  su  .-face 


H  +  X£  +  G  =  Q, 


(16) 


where  H  Is  the  sensible  heat  flux.  X£  the  latent  heat  flux,  and  G  the  9oll  heat  flux  which  Is  set 
equal  to  0.1  0-  For  snow-covered  surfaces.  G  Is  set  to  0  to  reflect  the  poor  conductivity  of  snow. 

DeBruln  and  Holtslag16  present  a  simplified  parameterization  for  the  sensible  heat  flux  (K) 
and  latent  heat  flux  (E) 


H  «  ■1.-Ta  +  r/l(Q-G)-(VY;  (17) 

l  +  y/s 

E  a  — £&—  (Q  -  G)  +  p'  y.  (18) 

1  +  y/s 


Here  s  *  q^/  T.  where  q^  Is  the  saturation  specific  humidity,  y  =  Cp/X  where  Cp  Is  the  specific 
heat  of  air  at  constant  pressure  and  X  the  latent  heat  of  water  vaporization,  and  a  and  P'  are 
empirical  parameters.  Holtslag  and  Van  Ulden  present  y/s  valueo  as  a  function  of  temperature 
for  a  standard  pressure  of  1000  mb.  A  regression  line  can  be  fitted  to  the  (^ata  such  that 


y/s  o  1 19.56  -  0.7843T  +  1 .288?  x  10-3  T2, 


(19) 


10  DeBruln,  H.A.R.  and  Holtslag,  A.A.M.  (1982)  A  simple  parameterization  of  the  surface 
fluxes  of  sensible  and  latent  heat  during  daytime  compared  with  the  Henman- Montelth 
concept,  J.  AppL  Meteorol.  21:1610-1621. 
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where  T  Is  the  air  temperature  In  Kelvin. 

From  Holtslag  and  Van  Ulden, 14  the  surface  moisture  parameter  a  Is  set  at  1 .0  for  wet 
conditions  and  0.45  for  dry  conditions.  The  parameter  p'  is  set  at  20  W  m'2. 

For  nighttime,  the  sensible  heat  flux  Is  related  to  the  total  cloud  cover  (N)  by  Smith  s17 
formula: 


H»  -40  (1  -N) 
where  H  Is  In  W  m2. 


(20) 


2.2.3  CALCULATING  SOLAR  ELEVATION  ANGLE 

The  solar  elevation  angle  for  a  given  time  and  location  may  be  calculated  by  a  method 
described  by  Woolf.18  The  following  relation  is  used  to  determine  the  solar  elevation  angle.  4>: 


sin  $  >  sin  LA  sin  D  +  cos  LA  cos  D  cos  H 


(21) 


where  LA  is  the  station  latitude,  D  the  solar  declination  angle,  and  H  the  solar  hour  angle. 

The  solar  declination  angle  Is  a  sinusoidal  function  of  time  with  maximum  and  minimum 
angles  occurring  during  summer  and  winter,  respectively.  There  is,  however,  a  slight 
asymmetry,  due  to  the  elllptlclty  of  the  earth's  orbit,  which  is  accounted  for  in  the  following 
expression  for  calculating  the  declination: 

sin  D  ■  sin  23.4438  sin  o,  (22) 


where 


o  (deg)  ■  a  +  279.9348  +  1.914827  sin  a  -  0.079525  cos  a 


(23) 


+  0.019938  sin  2a  -  0.001620  cos  2a. 


17  Smith,  F.B.  (1972)  A  scheme  for  estimating  the  vertical  dispersion  of  a  plume  from  a 
source  near  ground  level,  Proceedings  of  the  Third  Meeting  of  the  Expert  Panel  on  Air  Pollution 
Modeling.  NATO  committee  on  the  Challenges  of  Modem  Society,  Parts,  France. 

18  Woolf,  H.M.  (1980)  On  the  Computation  of  Solar  Elevation  Angles  and  the  Determln 

atlon  of  Sunrise  and  Sunset  Times,  National  Meteorological  Center,  Environmental  Sciences 
Services  Administration,  Hillcrest  Heights,  MO. 
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The  angular  fraction  of  a  year  (a)  for  a  particular  date  Is  given  by 


a  *  360  (JO-lJ/365.242, 


(24) 


where  JO  is  the  Julian  date. 

The  solar  hour  angle  (H).  a  measure  of  the  longitudinal  distance  to  the  sun  from  the  point 
for  which  the  calculation  Is  made,  is  given  by 


H(deg)  =  15(ZO  -  M)  -  LO 


(25) 


where  ZO  (hour)  Is  the  Greenwich  mean  time  (GMT)  of  the  calculation.  M  (hour)  the  time  of 
meridian  passage,  or  true  solar  noon,  and  LO  (deg)  the  station  longitude,  positive  being  west  of 
Greenwich.  M  Is  derived  from 


M  =  12  +  0.12357  sin  a  -  0.004289  C03  a  +  0.153809  sin  2a 
+  0.06078  cos  2a. 


(26) 


2.2.4  CALCULATING  THE  STABILITY  PARAMETER 
Method  1 

Based  on  data  collected  at  flve  sites,  Golder9  developed  a  nomogram  relating  the  Monln- 
Obukhov  length  (L)  and  surface  roughness  (ZJ  to  the  Fasquill  stability  categories.  This 

nomogram  Is  shown  as  Figure  1.  To  use  this  nomogram  In  the  model.  It  was  necessary  to  find 
a  mathematical  relationship  between  L,  ZQ,  and  the  stability  parameter  (SP).  The  following 
formulation  gives  a  good  approximation  of  SP  for  a  given  L  and  ZQ,  where  ZQ  is  the  surface 
roughness  at  the  wind  measurement  site. 


SP  =  A  +  B  loglo(100  ZJ 


(27) 


where  A  =  3.5  +  21.67/L 

B  a  0.48  when  11/Ll  >  0.015 
B  =  43.6311  /LI1  08  when  1 1/L|  <0.015 


B  ■  -B  when  1/L  <  0 
where  Z0  and  L  are  In  m. 


To  stay  within  the  bounds  of  the  chart,  Z0  is  restricted  to  the  0.1  to  40  cm  range  in  the 
model.  Values  of  SP  are  allowed  to  vary  from  0.5  to  6.0.  Extrapolating  below  0.5  is  not  very 
reliable,  and  using  0.5  as  a  minimum  value  is  a  conservative  approach. 

When  the  stability  parameter  is  equal  to  6.  which  occurs  at  nighttime  under  light  wind 
conditions,  the  model  assumes  that  horizontal  meandering  occurs.  When  concentration 
averaging  times  are  greater  than  1  min.  the  model  assumes  a  stability  parameter  for 
horizontal  diffusion  of  3.5,  or  a  neutral  condition.  Mean  concentrations  taken  over  averaging 
times  of  1  min  or  less  are  considered  to  be  unaffected  by  the  rather  large  oscillation  period  of 
the  meandering  wind.  Since  meandering  does  not  affect  the  vertical  diffusion,  the  stability 
parameter  for  the  vertical  diffusion  remains  at  6. 

If  the  lnverslc  >  V  3  less  than  50  m  above  the  ground  during  either  daytime  or 
nighttime,  the  s*«; .  .  parameter  for  the  vertical  diffusion  is  set  at  6.0.  At  night,  when 
meandering  ro&.  be  occurring,  the  stability  parameter  for  horizontal  diffusion  is  set  at  3.5 
when  the  concentration  averaging  time  is  >  1  min. 


Method  2 

If  the  standard  deviation  of  the  horizontal  wind  direction  (oe)  is  known,  then  the  stability 

parameter  can  be  calculated  using  the  Modified  Sigma  Theta  (MST)  approach  of  Mitchell.1 1  In 
this  method,  the  stability  categoiy  relationship  outlined  In  the  NRC  Regulatory  Guide  1.2319  is 
used.  However,  for  nighttime  when  the  ofl  is  greater  than  12.5  degrees.  Mitchell  assumes  that 

the  large  o#  is  due  to  meandering  and  not  to  Instability.  Since  meandering  does  not  affect  the 

vertical  diffusion,  an  adjustment  is  made  in  the  stability  classes  for  the  vertical  diffusion. 

The  MST  method  is  summarized  in  Table  3. 


19  USNRC  (1972)  Onsite  Meteorological  Programs  Regulatory  Guide  1.23,  U.S.  Nuclear 
Regulatory  Commission. 
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Figure  1.  1/L  as  a  Function  of  Pa  squill  Stability  Classes  and  Surface  Roughness  (zo). 


Table  3.  Modified  Sigma  Theta  (MST)  Method  to  Determine  Atmospheric  Stability  Class 
(from  Mitchell11) 


(Degrees) 

Daytime 

Stability 

Class2 

Wind  Speed 
(m  s'1) 

Nighttime1 

Stability 

Class3 

a6  >  22.5 

A 

u  <  2.4 

G 

2.4  <  u  <  2.9 

F 

2.9  <  u  <  3.6 

E 

3.6  <u 

D 

22.5  >  <Je  >  17.5 

B 

u  <2.4 

F 

2.4  <  u  <  3.0 

E 

3.0  <  u 

D 

17.5  >  ae  >  12.5 

C 

u  <  2.4 

E 

2.4  <u 

D 

12.5  >oe>  7.5 

D 

all  wlndspeeds 

D 

7.5>  oe  >  3.8 

E 

all  wlndspeeds 

E 

3.8  >  o0>  2.1 

F 

all  wlndspeeds 

F 

21  >«e 

G 

all  wlndspeeds 

G 

1  Nighttime  Is  defined  as  the  period  from  1  hour  before  sunset  to  1  hour  after  sunrise. 

2  More  applicable  to  describing  horizontal  dispersion  parameter  o  at  night. 

3  More  applicable  to  describing  vertical  dispersion  parameter  az  at  night. 


Aa  In  Method  1,  when  the  concentration  averaging  periods  are  1  min  or  less,  the 
concentrations  are  less  affected  by  the  rather  large  oscillation  period  of  the  meandering  wind. 
Therefore,  the  stability  parameter  for  the  horizontal  diffusion  is  set  equal  to  the  more  stable 
stability  parameter  for  vertical  diffusion. 

The  oe  values  In  Table  3  are  for  a  10  m  level  and  a  60-mln  time  period.  Therefore,  the 
measured  oe.  if  for  a  height  other  than  10  m  and  a  time  period  other  than  1  hour,  must  be 
adjusted  to  a  10  m  height  and  a  60-mln  time  period.  This  is  Qrst  done  by  adjusting  the  o0  to 
the  10  m  height  through  the  relationship 


(Oe)io-o,(IO/Z)02. 


(28) 


and  then  adjusting  the  10  m  oe  to  the  60-mln  time  period  using 
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(°»)eo  -  (o»)io  (60/t)0^. 


<29) 


where  Z  Is  the  wind  measurement  height  and  t  Is  the  time  period  over  which  oe  Is  measured. 

In  AFTOX,  the  following  equation,  derived  from  the  Information  In  Table  3.  is  used  to 
determine  the  stability  parameter  (SP)  for  all  cases  except  for  the  stability  parameter  for 
vertical  dispersion  at  night  when  c9  <  12.5°: 

SP  «  6.46  -  0.34 1  o8  +  0.0045  <J92.  (30) 


where  oe  Is  now  the  corrected  c9  for  a  10  m  height  and  60- min  time  period. 

Table  4  shows  how  the  model  determines  the  stability  parameters  for  vertical  diffusion  for 
nighttime  when  cb  •  12.5  degrees. 


Table  4.  The  Nighttime  Stability  Parameter  for  Vertical  Dispersion  as  a  Function  of  o9 
and  Wind  Speed 


(Degrees) 

Wind  Speed 
(m  sec'1) 

Vertical  Stability 
Parameter  (SP) 

o9>  22.5 

u  >  4.1 

SP  «  3.5 

2.4  <  u  <4.1 

SP  =  I4.4/u 

u  <  2.4 

SP  =  6.0 

22.5  >  c9  >  17.5 

u  >  3.4 

SP  =  3.5 

2.0  <  u  <  3.4 

SP  =  12/u 

u  <  2.0 

SP  =  6.0 

17. 5>  o0 

u  >2.7 

SP  ■=  3.5 

1.9  <u  <2.7 

SP  =9  ,6/u 

u  <  1.9 

SP  =  5.0 

The  AFTOX  model  makes  three  exceptions  to  the  MST  method: 

1)  During  the  daytime,  no  matter  how  small  the  oes  are.  the  model  does  not  allow,  as  In 
Method  1.  the  stability  p,irameters  to  exceed  3.5  (D  stability  class).  In  other  words,  only 
unstable  or  neutral  conditions  exist  during  the  daytime. 

2)  Any  time,  day  or  night,  that  the  height  of  the  inversion  base  Is  less  than  50  meters 
above  ground,  the  stability  parameter  for  the  vertical  diffusion  is  adjusted  to  6.0. 

3)  As  In  Method  1,  during  meandering  wind  or  low  Inversion  conditions,  the  stability 
parameter  for  the  horizontal  diffusion  is  set  equal  to  the  more  stable  stability  parameter  for 
vertical  diffusion  for  concentration  averaging  periods  of  1  min  or  less. 
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3. 3  Dispersion  Parameter* 

The  PaaqulU-Oifford  (P-G)  dispersion  parameters8  (oy.  o^)  are  used  m  the  dispersion  model. 
These  values  were  derived  from  the  Prairie  Grass  experiments 1,2-3  and  therefore  are 
appropriate  for  a  surface  roughness  of  3  cm,  and  are  based  on  10-mln  concentration  averages. 
The  oy  and  az  parameters  are  adjusted  for  other  surface  roughnesses  at  the  spill  site  (Z2q)  and 
concentration  averaging  times  (tjJ.  A  schematic  showing  the  steps  taken  In  the  model  to  arrive 
at  the  appropriate  dispersion  parameters  is  presented  In  Figure  2. 


Figure  2.  Procedure  for  Computing  ay  and  az 


2.3.1  HORIZONTAL  DISPERSION  PARAMETER  (o^ 

The  (P-G)  horizontal  dispersion  parameter  for  different  stabilities  and  for  any  distance  (x) 
can  be  approximated  by  using  a  power  law  curve 

oy*axb  (31) 
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where  a  and  b  are  given  In  Table  5  (from  Hansen2®  )  for  the  different  stability  categories  and 
parameters.  These  values  apply  when  oy  and  x  are  in  meters. 


Table  5.  Coefficients  and  Exponents  for  the  Horizontal  Dispersion  Parameter 


Stab.  Category(SC) 

A 

B 

C 

D 

E 

F 

Stab.  Parameter(SP) 

0.5 

1.5 

2.5 

3.5 

4.5 

5.5 

a 

0.40 

0.32 

0.22 

0.143 

0.102 

0.076 

b 

0.9 

0.9 

0.9 

0.9 

0.9 

0,9 

The  relationship  between  the  coefficient  a  and  the  stability  parameter  (SP)  can  be 
approximated  by  the  formula: 


a  =  0.479  -  0.1232  SP  +  0.00904  SP2.  (32) 

Since  Eq.  (31)  gives  oy  for  a  surface  roughness  of  3  cm  and  a  concentration  averaging  time 
of  10  min.  the  oy  must  be  corrected  for  other  surface  roughnesses  at  the  spill  site  (ZZJ  and 
pollutant  concentration  averaging  time  (tc).  The  following  relationship  Is  used  In  the  model  to 
arrive  at  the  corrected  horizontal  dispersion  parameter,  (oy)e: 


(°y)c  *  oy  (ZZ0/3)°  2  (t,  '10)°  2. 


(33) 


2.3.2  VERTICAL  DISPERSION  PARAMETER  (oz) 

The  vertical  dispersion  parameter.  oz,  for  a  roughness  length  of  3  cm  can  be  determined 
using  the  power  law 


o*  s  c  Jt*.  (34) 

where  x  and  oz  are  in  meters  and  c  and  d  are  given  In  Table  6.  The  constants  c  and  d  were 
derlvea  by  fitting  power  law  curves  to  the  Pasqulli  oz  curves. 


20  Hansen,  F.V.  (1979)  Engineering  Estimates  for  the  Calculation  of  Atmospheric 
Dispersion  Coefficients.  Internal  Report.  U.S.  Army  Atmospheric  Science  Laboratory.  White 
Sands  Missile  Range,  NM. 
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TABLE  6.  Coefficients  and  Exponents  for  the  Vertical  Dispersion  Parameter 


Stab.  Category(SC) 

A 

B 

C 

D 

E 

G 

Stab.  Parameter(SP) 

0.5 

1.5 

2.5 

3.5 

4.5 

5.5 

x(m)< 

745 

745 

2000 

1100 

1400 

1400 

c 

0.0414 

0.1036 

0.1173 

0.0975 

0.1050 

0.0617 

d 

1.3155 

1.0026 

0.9112 

0.8414 

0.7692 

0.7884 

x(m)> 

745 

745 

2000 

1100 

1400 

1400 

c 

1.928x1 O'4 

0.0534 

0.4422 

0.6097 

0.8788 

0.9990 

d 

2.1234 

0.7382 

0.5808 

0.4771 

0.4771 

The  oE  for  the  derived  stability  parameter  is  determined  by  linearly  interpolating  between  the 
ozs  for  the  stability  parameter  in  Table  6  Just  above  and  below  the  derived  stability  parameter. 
This  approach  was  taken  because  a  reasonable  line  of  best  fit  could  not  be  found  for  c. 

Through  the  numerical  solution  of  the  two-dimensional  equation  of  diffusion,  F.B. 

Smith17  derived  correction  factor  curves  for  oz  for  different  roughness  lengths.  Quite  often,  a 

simple  empirical  approximation  to  Smith's  results  is  used. 


Wc 


ot(ZZo/3\ 


,0.2 


(35) 


When  using  oe  to  determine  the  stability  parameter  (Method  2),  az  needs  to  be  corrected 
only  if  the  surface  roughness  at  the  spill  site  (ZZ^  is  different  than  the  surface  roughness  at 
the  wind  measurement  site  (ZJ.  The  corrected  vertical  dispersion  parameter.  (oz)c  is 


(0  )  =  o  (ZZ  /Z  ] 

'  z'c  z  '  o'  a 


,0.2 


(36) 


Correction  of  oz  for  the  concentration  averaging  time  is  not  needed  for  either  method  since 

large  eddies  are  much  less  likely  to  occur  in  the  vertical  direction  than  in  the  horizontal.  In 
the  model,  Eqs  (33),  (35)  and  (36)  apply  for  values  of  ZZQ  between  0.5  and  100  :m.  The  accuracy 

of  these  equations  is  questionable  beyond  these  limits. 
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2.4  Evaporation 


During  a  liquid  release,  it  Is  assumed  that  a  liquid  pool  will  be  formed  on  the  ground  after 
the  chemical  exits  the  original  container.  AFTOX  4.0  contains  three  evaporation  models.  The 
model  that  Is  used  depends  on  the  available  chemical  data.  The  three  models  are  1)  the  Vossler 
steady-state  energy  balance  model:21  2)  the  Shell  Oil  Company's  evaporation  model  used  in 
their  SPILLS  model;22  and  3)  the  Clewell  evaporation  model.23  For  a  more  detailed  description 
of  these  models,  the  reader  should  refer  to  the  above  references. 

2.4.1  VOSSLER  EVAPORATION  MODEL 

This  model  19  a  steady-state  model  in  which  the  pool  temperature  is  adjusted  based  on  the 
net  energy  input  into  the  pool  from  all  possible  sources.  The  pool  temperature  19  used  to 
determine  the  evaporation  rate.  These  sources  of  energy  include  9olar  radiation  (Qsol).  long 
wave  radiation  emitted  by  the  pool  (Q^l  and  the  atmosphere  (Qatm).  convective  heat  transfer 
from  the  atmosphere  (Qh<J,  heat  conducted  from  the  ground  (Q^),  and  heat  loss  due  to 
evaporation  (Q^).  The  steady-state  temperature  is  the  temperature  at  which  the  sum  of  all 
sources  of  heat  (energy)  transported  Into  the  pool  exactly  balances  the  heat  transfer  out  of  the 
pool:  that  is.  the  sum  of  energy  terms  described  above  is  zero.  The  9teady-state  energy  balance 
Is  expressed  as: 


Q«ol  +  Qatm  +  Qpol  +  Qhc  +  Qev  +  Qgrd  =  9  total' 


(37) 


where  Qtotaj  *  0  at  steady  state.  Many  of  these  energy  terms  can  be  expressed  as  a  function  of 
the  pool  surface  temperature.  The  equation  is  solved  iteratively  for  the  pool  temperature, 
which  is  then  used  to  calculate  evaporation  rate.  The  evaporation  rate  is  proportional  to  the 
mass  transfer  coefficient  at  the  liquid  pool-atmosphere  Interface  and  the  vapor  pressure  of  the 
chemical,  both  of  which  are  functions  of  the  pool  surface  temperature. 

This  model  requires  a  great  deal  of  chemical  data  as  can  be  seen  in  Section  3.1.2.  In  the 
model,  these  data  are  stored  In  a  separate  file  called  EVAP.DAT.  Currently,  only  four 
chemicals  are  stored  in  this  file.  They  include  nitrogen  tetroxide,  hydrazine, 
monomethylhydrazlne  (MMH).  and  u  ,mmetrical  dimethylhydrazine  (UDMH). 


21  Vossler.  T.L.  (1989)  Comparison  o/  Steady  State  Evaporation  Models  Jor  Toxic  Chemical 
Spills:  Development  of  a  New  Evaporation  Model  GL-TR-89-0319,  ADA221752. 

22  Fleischer,  M.T.  (1980)  SPILLS •  An  Evaporation/ Air  Dispersion  Model  for  Chemical 
Spills  on  Land.  Shell  Development  Company,  PB  83109470. 

23  Clewell,  H.J.  (1983)  A  Simple  Formula  for  Estimating  Source  Strengths  from  Spills  of 
Toxic  Liquids.  ESL-T-83-03. 
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2.4.2  SHELL  EVAPORATION  MODEL 

In  the  Shell  model,  the  evaporation  rate  Is  primarily  a  function  of  the  forced  convection 
over  the  liquid  pool  due  to  the  wind.  It  Is  assumed  that  there  Is  no  heat  transfer  due  to 
evaporative  cooling  or  radiation,  and  that  the  pool  temperature  Is  the  same  as  the  air 
temperature.  The  evaporation  rate  (Q)  Is  equal  to: 


0  =  (Sh  D/L)  L2  pv. 


(38) 


where  Sh  is  the  Sherwood  number.  D  Is  the  molecular  diffusion  coefficient.  pv  Is  the  vapor 
density,  and  L  Is  the  pool  width.  For  turbulent  flows  (Re  >  320000),  the  Sherwood  number  is 
equal  to: 


Sh  =-  0.037  Scl/3  (Re0  8  -  15200], 


(39) 


For  laminar  flows  (Re  <  320000).  the  Sherwood  number  is 


Sh=- 0.664  Sc 1/3  Re 1/2 


(40) 


where  Sc  la  the  Schmidt  number  and  Re  Is  the  Reynolds  number.  The  Schmidt  number  is 
given  by  v/D.  and  the  Reynolds  number  by  uL/v.  where  v  Is  the  kinematic  viscosity  of  air  In 
m2/sec.  and  u  Is  the  wind  speed  at  10  m  above  ground  level  In  m/sec. 

The  molecular  diffusion  coefflcent,  D,  is  equal  to: 

3/2f(M4  +  M  )/  M4  Ml 1/2 

D-cT  '  'lA - - A — EbL -  (4i) 

P°ABr‘0 

where  the  binary  dlffuslvlty  (D)  is  in  cm2/sec,  T  *s  the  temperature  in  Kelvin.  MA  and  MB  are 
molecular  weights  of  the  chemical  and  of  air,  P  is  the  pressure  in  atmospheres,  o^  is  the  hard 
core  diameter  of  the  binary  pair,  and  £2D  is  the  collision  Integral  for  diffusion. 

The  collision  Integral  £1D  Is  calculated  from  the  expression 

flD  =  A,  (T*)  Bi  +  Aj  exp(-3jT*),  (42) 
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where  the  A,  and  B,  are  constants  obtained  by  curve  fitting  calculated  values  of  ftD  versus  T*. 
T*  is  the  temperature  reduced  by  the  potential  well  depth  T*  may  be  expressed 


T*=T/Eab,  (43) 

where 

eAB  “  ^eA  ^2'  (44) 

where  e0  for  air  is  =  78.6.  The  AFTOX  chemical  data  file  for  the  chemical  Includes  eA  which  is 
approximately  equal  to  1.15  times  the  boiling  point  of  the  chemical. 

The  constants  in  Eq.  (42)  used  in  the  model  are  given  in  Table  7. 

Table  7.  The  A  and  B  Constants  Used  in  Eq.  (42)  in  AFTOX 


In  Eq.  (41)  oAB  is  equal  to 

°ab  =  0  5  <°a  +  °B>-  (45) 

where  oA  is  the  effective  diameter  of  the  gas  molecule  and  o0  Is  the  effective  diameter  of  the  air 
molecule  (set  at  3.71 1  in  the  model).  oA  is  included  in  the  AFTOX  chemical  data  file. 

The  constant  c  in  Eq.  (41)  is  taken  to  be 

103c  =  2,17  -  0.5|(Ma  +  Mb)/MaM011/2.  (46) 

The  molecu'ar  weight  of  air.  M0,  is  set  a»  28.966, 
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2.4.3  CLEWELL  EVAPORATION  MODEL 


For  those  chemicals  for  which  only  the  molecular  weight  and  vapor  pressure  are  known, 
the  Clewell  evaporation  model  Is  used.  This  model  is  a  crude  approximation  to  the  Ule  and 
Springer  model,24  which  is  similar  to  the  Vossler  model.  The  model  Is  represented  by  a  single 
formula: 


Q  «  0.08  u3/4  A  Ft  Z. 


(47) 


where  0  is  the  evaporation  rate  In  kg/hr,  u  Is  the  wind  speed  In  m/sec.  A  is  the  spill  area  In 
m2.  Ft  Is  a  pool  temperature  factor,  and  Z  Is  a  volatility  factor.  The  pool  temperature  factor  Is 

equal  to: 


F,  =  1  +  4.3  x  10'3  Tp2  for  Tp  >  0°C,  or  -  1  for  Tp  <  0°C.  (48) 

where  Tp  Is  the  pool  temperature  In  °C.  An  estimate  of  the  pool  temperature  must  be  made.  In 
general,  under  sunny  conditions  the  pool  temperature  may  be  several  degrees  warmer  than  the 
air  temperature,  and  for  highly  volatile  chemicals  It  may  be  several  degrees  cooler  than  the  air 
temperature. 

The  volatility  factor  (Z)  is  the  ratio  of  the  vapor  pressure  (VP)  and  molecular  weight  (MW) 
of  the  chemical  to  that  of  hydrazine,  or 


Z  *  (VP  MW)/(VPh  MWh)  . 


(49) 


In  the  model,  the  molecular  weight  for  hydrazine  Is  set  at  32.045  and  vapor  pressure  for 
hydrazine  is  set  at  .02113  atm. 


2.5  Saturated  Vapor  Pressure 

All  three  methods  of  determining  evaporation  rates  require  the  saturated  vapor  pressure 
(VP).  Depending  on  the  available  chemical  data,  one  of  two  formulae  is  used  to  compute  the 


34  Ille,  G.  and  Springer.  C.  (1978)  The  Evaporation  and  Dispersion  of  Hydrazine 
Propellants  from  Ground  Spills,  CEEDO-TR-83-30,  ADA  059407. 
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vapoi  pressure.  They  are  the  Antoine  equation  and  the  Frost- Kalkwarf  equation.  The  Antoine 
equation  is: 


VP-At-Br/fTp  +  CfJ. 


(50) 


where  A^,  Br.  and  Cr  are  empirical  constants,  and  are  Included  in  the  AFTOX  chemical  data 
file. 

The  FTost-Kalktvarf  equation  is: 


In  Pr  -  Ar(l  -  1/Tr)  -  Bf  In  Tf  «•  Dr(Pf/Tr  -  1). 


(51) 


where  Pf  is  the  reduced  vapor  pressure  VP/Pc,  Tf  is  the  reduced  temperature.  T  /Tc.  and  where 
Pc  and  Tc  are  the  critical  pressure  and  temperature,  respectively,  and  are  Included  in  the 
chemical  data  file.  Dr  is  equal  to  27/64,  and  A^.  and  Br  are  empirical  constants  Included  in  the 
chemical  data  Ole. 

The  vapor  pressure  is  then  equal  to: 


VP  -  Pr  Pc. 


(52) 


Either  equation  Is  used  for  determining  the  vapor  pressure  for  the  Shell  evaporation 
model.  Since  the  Shell  model  assumes  that  the  pool  temperature  is  equal  to  the  air 
temperature,  Tp  Is  the  air  temperature.  The  VossleT  model  uses  an  expanded  version  of  the 

Antoine  equation.  The  Clewell  model  uses  a  constant  vapor  pressure  that  is  either  stored  in 
the  chemical  data  file  or  entered  by  the  user  when  running  the  model  for  a  chemical  not  in  the 
chemical  data  file. 


2.0  Liquid  Density 

When  the  source  strength  of  a  liquid  release  is  given  in  volume  per  unit  time  for  a 
continuous  spill  or  volume  for  an  instantaneous  spill,  the  liquid  density  must  be  known  to 
convert  the  source  strength  to  mass. 

The  liquid  density  can  be  calculated  from  the  Guggenheim  equation 


pL-MW|l  +a(l  -Tr),/3  +  b  (1  -Tr)]/Vc, 


(53) 
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where  a  and  b  are  empirical  constants  that  art  Included  in  the  chemical  data  file.  When  the 
constants  are  not  available,  a  liquid  density  value  may  be  Included  in  the  file  In  place  of  the 
"a"  constant.  If  no  values  are  available,  the  model  assumes  a  liquid  density  of  1  g/cm3.  Vc  is 

the  critical  volume  and  is  Included  In  the  chemical  data  file. 


2.7  Continuous  Buoyant  Plumes 

The  buoyant  plume  routine  considers  a  steady  state  emission  from  a  single  stack  over  a 
finite  or  Infinite  time  period.  The  final  centerline  rise  height  and  the  distance  downwind  at 
which  the  plume  reaches  that  height  are  calculated  using  the  widely  accepted  Briggs  plume  rise 
formulas.26 

The  initial  buoyancy  flux  (F)  is  defined  by 

F-lg/n)  V(T.-Tj/TB.  (54) 


where  g  is  the  acceleration  of  gravity  (m/sec2).  T#  is  the  temperature  of  the  plume  (°K).  T  the 
ambient  temperature  (°K(.  and  V  the  plume  volume  flux  (ra3/sec). 

For  a  buoyant  bent-over  plume  In  unstable  or  neutral  conditions,  the  height  of  the 
centerline  of  the  plume  is  given  by 


Z«  1.6  F 1/3  X?/3/u 


(55) 


where  Z  and  X  ate  in  km  and  u  is  wind  speed  in  ra/sec. 

The  plume  reaches  Its  equilibrium  height  at  a  distance  (X*)  downstream  defined  as  the 
point  at  which  atmospheric  turbulence  dominates  the  entrainment.  In  the  model  this  distance 
is  defined  as  follows: 

For  F  >  55 

X*  -  1 19  F04;  (56a) 


26  Briggs,  G.A.  (1975)  Plume  Rise  Predictions,  Lectures  on  Air  Pollution  and  Environ¬ 
mental  Impact  Analyses,  Workshop  Proceedings,  Boston.  MA  Sept.  27  -  Oct.  3,  1975,  pp.  59-11 1, 
American  Meteorological  Society,  Boston,  MA. 
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For  F  <  55 


X*  -  49  F0-625. 


(56b) 


For  a  buoyant,  bent-over  plume  In  stable  air,  the  final  plume  centerline  height  is  given  by 


Z  =»  2.6 


(57) 


For  a  buoyant,  vertical  plume  with  a  wind  of  1  m/sec  or  less,  the  final  centerline  height  Is 
given  by 


Z  =  5.0  F1/4  s'3/®, 


(58) 


where  s  Is  a  stability  parameter  equal  to: 


s  =  g/T  DT/DZ. 


(59) 


and  where  DT/DZ  Is  the  potential  temperature  lapse  rate,  which  in  the  model  Is  set  equal  to 
0.02°/m  for  slightly  stable  and  0.035°/m  for  more  stable  conditions. 

The  distance  downwind  where  the  plume  reaches  the  equilibrium  height  Is  defined  as: 


X'-rtu/s1/2.  (60) 

Unlike  the  unstable  and  neutral  case,  this  distance  is  not  required  for  computing  the 
equilibrium  height,  but  is  calculated  and  printed  out  for  information  purposes  only. 

The  source  effective  emission  height.  H,  Is  equal  to  the  sum  of  the  equilibrium  height.  Z, 
and  the  stack  height.  Hs: 


H  -  H,  +  Z. 


(61) 
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If  H  is  higher  than  the  inversion  height,  then  H  is  set  equal  to  the  inversion  height.  In  other 
words,  we  assume  that  the  buoyant  plume  does  not  penetrate  the  inversion.  This  is  a 
conservative  approach  when  calculating  concentration  levels  near  the  ground.  In  calculating 
the  distribution  of  pollutant  concentrations,  the  model  assumes  that  the  pollutants  are 
released  at  the  source  effective  emission  height  and  not  at  the  stack  height. 


3.  AFTOX  COMPUTER  PROGRAM  OPERATION 


This  section  provides  the  program  structure  and  execution  instructions  for  the  AFTOX 
computer  program.  The  code  is  written  in  the  GW-BASIC  programming  language. 


3.1  File  Structure 

The  AFTOX  program  contains  five  component  program  files  and  six  data  Bles,  four  for 
input  data  and  two  for  output.  Three  ancillary  files  are  used  for  editing  the  data  in  three  of  the 
input  data  files.  CHAIN  statements  are  used  to  connect  the  various  program  81es.  A  list  of  the 
files  is  shown  in  Table  8. 
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Table  8.  List  of  files  In  AFTOX 


1  PROGRAM  FILES  | 

AFTOX.EXE 

Introduction 

DSP1.EXE 

Defines  chemical  properties  and  the  meteorological 
conditions 

DSP2.EXE 

Defines  source  conditions  (that  Is.  emission  rate, 
spill  duration,  and  spill  area) 

DSPHP.EXE 

Defines  source  conditions  for  buoyant  plume  from  a  | 
stack 

DSP3.EXE 

Computes  1)  the  hazard  area.  2)  the  maximum 
concentration,  and  3)  the  concentration  at  a  given  | 

point  and  time  | 

I  DATA  FILES  f 

SD.DAT 

Station  data  1 

CH.DAT 

Chemical  name  and  data,  and  toxic  limits 

EVAP.DAT 

Chemical  data  for  Vossler  evaporation  model 

AFT.  DAT 

File  for  storing  all  input  and  output  data  when 
printer  is  off 

CONCXY.DAT 

File  for  storing  x,  y  positions  of  concentration 
contour  | 

DEV1CE.DAT 

Contains  information  on  computer  and  screen  type 

1  ANCILLARY  FILES 

SETUP.EXE 

Establishes  computer  and  screen  type  | 

SDFIL.EXE 

Edits  station  data  file  (SD.DAT)  g 

CHFIL.EXE 

Edits  chemical  data  file  (CH.DAT)  1 

3.1.1  ANCILLARY  FILES 
SETUP.EXE 

Before  running  AFTOX  for  the  first  time,  SETUP.EXE  should  be  run  to  establish  the  type  of 
computer  and  screen.  When  running  SETUP.EXE,  a  screen  as  shown  In  Figure  3  will  appear. 
Enter  the  appropriate  code  for  the  computer  and  monitor  types.  The  appropriate  information 
will  be  stored  in  DEV1CE.DAT. 
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AFTOX 


AIR  FORCE  TOXIC  CHEMICAL  DISPERSION  MODEL 

SET  UP  FOR  GRAPHICS  PROGRAM 

Enter  code  for  type  of  graphics  card/monitor. 

E  ■  EGA/VGA  C  ■  CGA  N  ■  Other  or  no  graphics  ?  E 

Enter  code  for  monitor  type 
C  -  COLOR  M  ■  MONOCHROME  ?  C 


Figure  3.  Screen  display  for  SETUP  routine. 


SDFIL.EXE 

The  SDFIL.EXE  file  Is  used  to  set  up  the  station  data  file.  SD.DAT.  SDFIL  allows  one  to 
view,  delete,  add.  and  edit  data  in  SD.DAT.  The  Ole  may  contain  one  or  more  stations. 

To  run  SDFIL.  type  SDFIL  and  press  <EN1ER>.  The  computer  will  display  a  screen  with 
several  options  as  shown  in  Figure  4.  Choose  the  appropriate  option  and  proceed  as  Instructed. 


1—  display  or  print  station  data 

2- -set  up  new  SD.DAT  file 

3-  -edit  data 

4—  add  station 

5—  delete  station 

6- -qult 

Choose  one  of  the  above 


Figure  4.  List  of  options  In  SDFIL. 
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The  following  Is  a  list  of  data  required  for  each  station  stored  In  SD.DAT. 


Station  name 
1 -Metric  2-Engllsh  Units 
Standard  Deviation  of  Wind  Direction  (Y/N) 
Standard  Deviation  Averaging  Time  (min) 
Latitude  (deg) 

Longitude  (deg) 

Surface  Roughness  (cm) 

Height  of  Wind  Measurement  (m  or  ft) 

Time  Difference  (Greenwich-local  standard) 
Station  Elevation  (m  or  ft) 


Stations  that  are  located  reasonably  close  to  each  other  (<50  km)  do  not  have  to  be  listed 
separately  because  the  latitude,  longitude,  and  elevation  are  not  highly  sensitive  parameters. 
However,  if  the  stations  have  different  surface  roughnesses  or  wind  measurement  heights  then 
they  must  be  listed  separately. 

The  units  refer  to  the  meteorological  and  distance  parameters.  The  user  will  have  the 
option  of  entering  either  metric  or  English  units  for  the  spill  rate  or  quantity. 

As  a  guide  for  inputting  the  proper  surface  roughness,  a  table,  shown  below.  Is 
automatically  displayed  on  the  screen.  The  surface  roughness  need  not  be  restricted  to  one  of 
the  values  shown  in  Figure  5.  but  It  Is  recommended  that  the  value  fall  between  0.5  and  100 
cm. 


TERRAIN  DESCRIPTION _ SURFACE  ROUGHNESS  (cm) 


SNOW.  NO  VEGETATION.  MUD  FLATS.  NO  BUILDINGS  0.5 

RUNWAY.  OPEN  FLAT  TERRAIN.  GRASS.  FEW  SMALL  BUILDINGS  3 

LOW  CROPS.  OCCASIONAL  LARGE  BUILDINGS  10 

HIGH  CROPS.  SCATTERED  BUILDINGS  25 

PARKLAND.  BUSHES.  NUMEROUS  BUILDINGS  50 

REGULAR  LARGE  OBSTACLE  COVERAGE  (SUBURB,  FOREST)  100 


MM-tartr.  a 


Figure  5.  Roughness  length  as  a  function  of  tmalr^  type. 
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CHF1L.EXE 

The  CHFIL.EXE  file  is  used  to  view,  make  changes,  deletions,  or  additions  to  the  chemical 
data  file,  CH.DAT.  The  current  CH.DAT  file  contains  129  chemicals.  A  list  of  the  chemicals  Is 
given  in  the  Appendix.  The  chemicals  are  listed  In  alphabetical  order,  except  for  the  IJ.S.  Air 
Force's  CH.DAT  file  in  which  the  first  ten  chemicals  listed  are  chemicals  most  frequently  used 
by  the  Air  Force. 

To  run  CHF1L,  type  CHFIL  and  press  <ENTER>.  The  computer  will  display  a  screen  With 
several  options  as  shown  in  Figure  6.  Choose  the  appropriate  option  and  proceed  as  instructed. 


1 — display  or  print  data 

2—  edit  exposure  limits  for  specific  chemical 

3—  edit  data  for  specific  chemical 

4—  add  chemical 

15— delete  chemical 
6-quit 


Choose  one  of  the  above 


Figure  6.  list  of  options  In  CHFIL. 


In  viewing  the  chemical  data,  several  options  exist:  1)  data  may  be  displayed  on  the  screen 
or  sent  to  the  printer;  2)  either  the  chemical  numbers  and  names  or  all  the  chemical  data  may 
be  displayed  or  printed:  and  3)  chemical  data  for  a  particular  chemical  or  for  all  the  chemicals 
may  be  displayed  or  printed. 

The  exposure  limits  are  used  only  as  a  guide  for  the  user.  The  user  must  Input  the  exposure 
limits  when  running  AFTOX,  but  tnay  default  to  the  STEL  values. 

When  editing  the  chemical  data  under  option  3,  the  user  must  re-enter  all  the  data  for  a 
particular  chemical,  even  though  he  may  be  changing  only  one  value.  The  old  data  Is 
displayed  on  the  screen  to  make  it  easy  for  the  user  to  re-enter  the  data.  One  word  of  caution— 
the  names  hydrazine,  monomethylhydrazine  (MMH),  dlmethylhydraztne  (UDMH),  and 
nitrogen  tetxoxlde  must  be  spelled  exactly  because  the  model  must  recognize  the  name,  and  not 
the  number.  In  order  for  it  to  use  the  Vossler  evaporation  model. 

The  data  for  each  chemical  consist  of  the  following: 


Chemical  number 
Chemical  name 
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Time  weighted  average  (TWA)  exposure  limits  (PPM  and  mg/m3) 
Short  term  exposure  limits  (STEL)  (PPM  and  mg/m3) 

Molecular  weight 
Bolling  temperature  (K) 

Critical  temperature  (K) 

Critical  pressure  (atm) 

Critical  volume  (cm3/g-mole) 

Vapor  pressure  constant  (3) 

Liquid  density  constants  (2) 

Molecular  dlffuslvlty  constants  (2) 


If  the  chemical  is  a  gas  at  all  expected  ambient  temperatures,  then  only  information  up 
through  the  boiling  temperature  is  needed. 

If  the  chemical  is  a  liquid,  vapor  pressure  information  must  be  In  the  chemical  (lie.  The 
exception  is  for  the  four  chemicals,  hydrazine,  raonomethylhydrazlne  (MMH), 
dlmethylhydrazlne  (UDMH).  and  nitrogen  tetroxlde.  that  use  the  Vossler  evaporation  model. 

T1  ?  chemical  data  required  for  the  Vossler  model,  including  the  vapor  pressure  data,  are 
included  In  EVAP.DAT. 

The  vapor  pressure  may  be  defined  In  one  of  three  ways:  1)  by  the  Antoine  equation,  in 
which  case  three  constants  are  required,  2)  by  the  Frost- Kalkwarf  equation  in  which  two 
constants  are  required,  or  3)  by  directly  Inputting  the  vapor  pressure  In  atm.  The  model 
determines  which  method  to  use  by  the  number  of  constants  In  the  file.  If  the  Frost-Kalkwarf 
equation  is  to  be  used,  the  critical  temperature  and  critical  pressure  must  be  known. 

The  liquid  density  is  required  for  all  chemicals  In  the  liquid  state,  including  those  using 
the  Vossler  evaporation  model.  The  liquid  density  is  defined  in  one  of  two  ways:  1)  by  the 
Guggenheim  equation,  in  which  case  two  constants  are  required,  or  2)  by  directly  Inputting  the 
liquid  density  in  g/cm3.  As  with  the  vapor  pressure,  the  model  determines  which  method  of 
defining  the  liquid  density  by  the  number  of  constants.  The  Guggenheim  equation  requires 
both  the  critical  temperature  and  critical  volume.  If  no  density  is  entered,  the  model  will 
assume  a  density  of  1  g/cm3. 

The  molecular  dlffuslvlty  constants  are  the  effective  diameter  of  the  molecule  (A)  and  the 
energy  of  molecular  interaction  (J).  The  dlffuslvlty  is  required  for  the  Shell  evaporation 
model.  If  either  of  these  two  constants  is  not  available,  the  model  defaults  to  the  Clewell 
evaporation  model. 

3.1.2  DATA  FILES 

There  are  six  data  files,  three  of  which  have  been  discussed  in  Section  3.1.1.  The 
remaining  three  are  EVAP.DAT.  CONCXY.DAT.  and  AFT. DAT. 
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EVAP.DAT 


EVAP.DAT  is  the  chemical  data  file  for  the  Vossler  evaporation  model.  The  chemicals 
listed  In  this  file  are  hydrazine,  monomethylhydrazlne  (MMH),  dimethylhydrazlne  (UDMH), 
and  nitrogen  tetroxide.  The  chemical  data  Include: 


Molecular  weight 
Bolling  temperature 
Freezing  temperature 
Molecular  diffusion  volume 

Constants  and  applicable  temperature  ranges  for: 
Vapor  viscosity 
Vapor  heat  capacity 
Vapor  thermal  conductivity 
Heat  of  vaporization 
Saturation  vapor  pressure 
Liquid  thermal  conductivity 


Data  may  be  changed  using  an  edit  program.  However.  If  a  new  chemical  is  added  to  the 
file,  the  code  In  DSP2.BAS  must  be  changed  so  that  the  new  chemical  will  be  recognized. 

CONCXY.DAT 

CONCXY.DAT  Is  an  output  data  file  containing  Information  on  the  X.  Y  coordinates  of 
up  to  three  requested  concentrations.  Data  stored  In  the  file  include  the  concentration  of 
Interest,  either  mg/m3  or  ppm.  the  time  since  release  for  an  Instantaneous  or  finite  continuous 
release,  and  the  contour  half  width  for  various  distances  downwind  In  either  m  or  ft.  The  half 
width  has  20  m  resolution.  The  actual  half  width  will  be  within  20  m  but  always  less  than  the 
calculated. 


3.2  8ettlng  Up  AFTOX 

AFTOX  may  be  run  directly  from  the  diskette,  but  If  your  computer  has  a  hard  disk  it  Is 
recommended  that  the  program  be  transferred  to  the  hard  disk.  Create  a  directory  on  your 
hard  disk  where  you  would  like  the  program  to  reside  (for  example,  C:>md  AFTOX).  Then  copy 
the  flies  from  the  diskette  to  the  hard  disk  (for  example,  A:copy«.»  C:>AFTOX). 

To  tailor  the  program  to  your  particular  computer,  run  SETUP.EXE.  Specify  whether  your 
machine  Is  VGA/EGA,  CGA,  or  no  graphics.  If  your  machine  has  Hercules  graphics,  you  will 
have  to  choose  no  graphics.  Also,  specify  whether  your  monitor  is  color  or  monochrome. 
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Once  SETUP  has  been  run,  there  will  be  no  need  to  run  It  again  unless  you  transfer  the 
program  to  a  different  type  of  computer. 

Information  on  your  location  must  be  stored  In  the  station  data  file,  SD.DAT,  by  running 
SDFIL.EXE.  Once  you  input  the  station  data  into  SD.DAT.  there  will  be  no  need  to  run 
SDFIL.EXE  again  unless  you  change  your  location. 

You  are  now  ready  to  run  AJFTOX. 


3.3  Running  AFTOX 

3.3.1  GENERAL 

AFTOX  is  a  very  user-friendly  program.  Simply  proceed  through  the  program,  answering 
the  questions  as  you  go  along.  Default  values  are  frequently  given  <ln  brackets>.  If  you  wish 
to  go  with  the  default  value,  simply  press  <ENTER>.  The  program  has  the  feature  of  being  able 
to  back  up  by  entering  <999>  <ENTER>,  thus  eliminating  the  need  to  start  over  if  you 
accidently  entered  the  wrong  data.  The  exception  is  that  you  can  not  back  up  into  the  previous 
file  (for  example,  you  can  not  back  up  from  DSP2  Into  DSP1). 

3.3.2  PRINTING  AND  STORING  DATA 

Type  AFTOX  at  the  DOS  prompt  and  press  <ENTER>.  This  will  start  the  program  execution. 
If  you  are  using  a  printer,  and  wish  the  plume  plot  to  be  sent  to  the  printer,  you  must  load  the 
appropriate  PSC  utility  program  before  you  start  AFTOX.  If  you  haven't,  the  program  will 
remind  you.  and  you  will  have  to  start  over.  If  the  printer  Is  off,  then  the  Input  and  output 
data  are  stored  In  AFT.DAT  for  viewing  or  printing  out  at  a  later  time.  Each  time  AFTOX  Is 
run,  the  data  In  AFT.DAT  is  erased.  If  you  wish  to  save  the  data  In  AFT.DAT.  the  file  can  be 
saved  by  renaming  the  file  (for  example,  C:>AFTOX:ren  aft.dat  spill#  1  .dat). 


3.3.3  STATION  DATA 

If  there  is  more  than  one  station  In  the  station  data  file  (SD.DAT),  you  will  be  asked  to 
enter  the  appropriate  station.  The  data  stored  in  SD.DAT  has  been  discussed  In  Section  3.1.1. 

3.3.4  DATE  AND  TIME 

The  computer  date  and  time  are  displayed.  If  the  spill  is  for  a  different  date  or  time,  a  new 
date  and  time  may  be  entered.  The  program  converts  the  date  to  a  Julian  date  which  it  uses  In 
conjunction  with  the  time,  latitude  and  longitude  to  determine  the  solar  elevation,  and 
subsequently,  the  solar  insolation  and  surface  heat  flux. 
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3.3.5  TYPE  OF  RELEASE 


AFTOX  handles  five  types  of  releases: 


Continuous  gas 
Continuous  liquid 
Instantaneous  gas 
Instantaneous  liquid 
Continuous  buoyant  stack 


The  user  has  a  choice  of  a  continuous,  instantaneous,  or  buoyant  release.  An  Instantaneous 
release  Is  defined  as  occurring  over  a  15-sec  period.  A  continuous  release  Is  any  release 
occurring  o/er  a  period  greater  than  15-sec.  For  the  continuous  and  instantaneous  releases, 
the  model  determines  whether  It  is  a  gas  or  liquid,  based  on  whether  the  air  temperature  is 
above  or  below  the  boiling  point  of  the  chemical.  If  the  type  of  release  is  unknown,  the  user 
should  choose  an  instantaneous  release  for  a  conservative  corridor. 

3.3.6  CHEMICAL  DATA 

A  list  of  129  chemicals  is  displayed  on  the  screen.  This  fist  is  shown  in  Appendix  A.  The 
user  enters  the  appropriate  number.  If  the  chemical  of  interest  is  not  listed,  the  user  may  press 
<ENTER>  at  the  end  of  the  list  of  chemicals.  He  will  then  be  asked  the  name  of  the  chemical 
and  its  molecular  weight.  If  the  molecular  weight  is  entered,  the  model  asks  for  the  vapor 
pressure  in  mm  Hg.  If  the  vapor  pressure  is  known,  the  model  uses  Clewell's  lormula  for 
determining  the  evaporation  rate.  If  either  the  molecular  weight  or  vapor  pressure  ie  not 
known,  the  model  assumes  the  worst  case;  that  is,  the  evaporation  rate  is  equal  to  the  spill 
rate.  Also,  if  the  molecular  weight  is  not  entered,  the  concentrations  must  be  in  mg/m3  since 
conversion  to  ppm  is  not  possible  without  knowing  the  molecular  weight. 

For  a  buoyant  plume  release,  the  model  bypasses  the  chemical  list  and  data  file  and  asks 
only  for  the  molecular  weight.  Again,  if  the  molecular  weight  is  not  known,  concentrations 
must  be  in  mg/m3. 

3.3.7  METEOROLOGICAL  DATA 

The  meteorological  data  consist  of: 


Air  temperature 
Wind  direction 
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Wind  speed 

Standard  deviation  of  wind  direction  and  time  over  which  it  is  determined  (optional) 
Cloud  amount 
Predominant  cloud  category 

Ground  condition  -  iry,  wet,  snow  covered  (daytime  only) 

Inversion  base  helg.it 


The  air  temperature  Is  a  necessary  input  parar  tetei  but  does  not  have  a  large  Influence  on 
the  results.  If  a  temperature  reading  is  not  avallab  e.  a  reasonable  guess  would  be  sufficient. 

Calm  winds  are  not  allowed.  If  a  zero  wind  spee^’  lf>  entered,  the  model  will  adjust  the  speed 
to  either  0.5  m/sec  or  1  kt.  depending  on  which  units  th»  user  chooses.  The  model  converts  the 
wind  speed  measurement,  whose  height  Is  specified  In  the  station  file  SD.DAT,  to  a  10-ni 
height  wind  speed.  The  1  »-m  wind  speed  Is  used  In  all  of  the  calculations,  and  therefore,  the 
user  should  be  aware  that  the  plume  may  move  down  wind  at  a  faster  rate  than  the  measured 
wind  speed  wo**M  indicate. 

If  the  standard  deviation  of  wind  direction  Is  normally  available,  it  Is  to  Indicated  In 
SD.DAT  along  with  the  time  over  which  it  Is  determined.  If  the  standard  deviation  Is  not 
available  at  the  time  that  the  model  Is  being  run.  simply  press  <ENTER>  and  the  model  will 
default  to  using  the  wind  speed  and  solar  conditions  for  computing  stability  and  the  corridor 
width. 

The  cloud  amount  is  entei  ;d  in  elghuis.  There  are  three  cloud  categories  to  choose  from  - 
high,  middle,  and  low.  If  there  are  two  cloud  layers  present,  the  operators  should  use  the  layer 
with  the  larger  cloud  amount.  If  there  ar;  two  layers  with  equ£  !  cloud  amount,  the  operator 
shou.M  choose  the  lower  cloud  layer.  The  cloud  type  Is  not  a  factor  at  night. 

There  ere  three  ground  types  to  choose  from  -  wet,  dry.  snow  covered.  If  ir.  doubt  as  to 
whether  the  ground  Is  wet  or  dry,  the  user  should  choose  wet,  which  will  result  In  a  more 
conservative  hazard  distance.  If  the  air  temperature  Is  20°C  (68°F)  or  greater,  the  model 
assumes  no  snow  cover. 

If  the  base  of  ar.  inversion  Is  below  500  m  (2000  ft),  the  height  of  its  base  Is  entered. 
Inversion  heights  greater  than  500  m  (2000  ft)  may  be  entered  but  will  have  no  ef  ect  on  die 
plume.  If  the  base  of  the  Inversion  Is  below  50  m  (lf4  ft),  the  model  assumes  no  Inversion  but 
does  assume  a  stability  parameter  of  6  When  ther<;  is  an  inversion,  the  model  assumes  that 
the  pollutants  are  trapped  below  its  base.  In  the  rare  case  where  the  release  is  ahnve  the 
Inversion,  the  pollutants  remain  above  the  Inversion. 

3.3.8  ROUGHNESS  LENGTH  AT  SPILL  SITE 

For  a  non-buoyant  release,  the  surface  roughness  length  at  the  spill  site  must  be  entered. 

For  a  buoyant  plume.  It  Is  assumed  that  the  elevated  plume  is  minimally  affected  by  the 
surface  roughness,  which  is  set  at  3  cm.  The  user  has  the  option  to  call  up  the  table  of 
roughness  lengths  (Figure  5)  for  different  types  of  terrain.  When  uncertain  as  to  the 
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appropriate  roughness  length  for  the  spill  site,  the  user  should  choose  a  lower  value,  which 
will  produce  longer  hazard  distances^  Input  roughness  lengths  below  0.5  cm  or  greater  than 
100  cm  are  set  at  0.5  or  100  cm.  respectively. 

3.3.9  SOURCE  INFORMATION 

The  source  Information  required  varies  depending  on  the  type  of  spill  -  continuous. 
Instantaneous,  gas,  or  liquid. 

Continuous  Gas  Release 

For  a  continuous  gas  release,  the  following  information  Is  required: 


Height  of  leak  above  ground  (m.  ft) 

Emission  rate  through  rupture  (kg/mln,  lb/min) 
Total  time  of  release  (min) 


The  emission  rate  Is  assumed  constant  for  the  total  time  of  the  spill.  If  the  duration  of  the 
release  is  finite  then  the  total  amount  released  is  displayed. 

Continuous  Liquid  Release 

For  a  continuous  liquid  release,  the  following  information  is  required: 


Spill  rate  through  rupture  (kg/mln,  m3/mln,  lb/mln,  gal /min) 

Total  time  of  spill  (min) 

Spill  area  (m2,  ft2),  otherwise  dt  alue  19  used 

Pool  temperature,  for  those  chemlc  )  using  the  Clewell  evaporation  model  (default  «*  air 
temperature) 


The  height  of  the  leak  is  not  entered  since  it  is  assumed  that  the  liquid  spills  to  the  ground, 
"he  default  spill  area  Is  based  on  the  volume  spilled  and  the  assumption  of  a  1-cm  deep  pool. 
For  an  ongoing  continuous  spill,  the  volume  spilled  Is  based  on  a  10-mln  spill.  The  default 
spill  area  Is  displayed  but  if  the  user  has  information  on  the  size  of  the  spill  area,  he  may 
over-ride  the  default  value.  Based  on  the  Input  data  and  the  chemical  properties,  the  model 
computes  the  evaporation  (emission)  rate  Into  the  atmosphere  and  the  total  evaporation 
(release)  time  for  a  finite  continuous  release.  The  model  assumes  that  the  spill  area  remains 
constant  with  time. 
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Instantaneous  Gas  Release 


For  an  Instantaneous  gas  release,  the  following  Information  is  required: 


Release  height  (m,  ft) 
Amount  released  (kg,  lb) 


The  model  assumes  a  cylindrical  volume  source  in  which  the  height  is  equal  to  the  radius. 
The  Initial  volume  of  the  spill  is  a  function  of  the  amount  released. 

Instantaneous  Liquid  Release 

For  an  instantaneous  liquid  release,  the  following  information  is  required: 


Amount  spilled  (kg,  m3.  lb.  gal) 

Spill  area  (m2,  ft2),  otherwise  default  value  Is  used 

Pool  temperature,  for  those  chemicals  using  the  Clewell  evaporation  model  (default  ■  air 
temperature) 


The  default  spill  area  is  based  on  the  volume  spilled  and  the  assumption  of  a  1-cm  deep 
pool.  Based  on  the  input  data  and  the  chemical  properties,  the  model  computes  the  emission 
rate  into  the  atmosphere  and  the  total  time  of  release. 

3.3.10  WORST  CASE  SCENARIO 

When  the  initial  spill  alert  is  sounded,  quite  often  very  little  source  information  is 
available.  In  this  case.  AFTOX  4.0  has  the  option  of  computing  a  worst  case  scenario.  The 
toxic  corridor  length  calculation  Is  based  on  Air  Force  Regulation  355-1/AWSSUP1. 
Attachment  1.  dated  24  September  1990.  This  regulation  states  that  the  corridor  length  is 
equal  to  the  current  wind  speed  In  knots  x  6,000.  This  represents  the  distance  in  feet  that  the 
plume  will  travel  In  one  hour. 

The  corridor  width  Is  a  function  of  the  wind  speed,  or  standard  deviation  of  wind  direction 
if  available.  The  computaUon  of  the  width  Is  described  In  Section  3.3.15.  The  exception  is 
that  there  Is  no  width  adjustment  for  the  duration  of  the  spill  since  the  duration  may  be 
unknown. 

Having  computed  the  worst  case  sctnarlo,  the  user  can  then  proceed  and  enter  the  source 
Information  when  It  becomes  available  without  re-entering  the  meteorological  information. 
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However,  if  you  need  to  change  the  name  of  the  chemical  or  the  source  type  (continuous, 
instantaneous),  then  you  must  start  over 

3.3.11  CONTINUOUS  BUOYANT  PLUME  RELEASE 

For  a  continuous  buoyant  plume  release  horn  a  stack,  the  following  Information  is 
required: 


Molecular  weight  (if  available) 
Emission  rate  (kg/min  or  lbs/mln) 
Elapsed  time  of  emissions  (min) 

Stack  height  (m  or  ft) 

Gas  stack  temperature  (°C  or  °F) 
Volume  flow  rate  (m3/min  or  ft3/min) 


With  these  data,  the  model  calculates  the  equilibrium  plume  height  and  the  downwind 
distance  at  which  the  plume  reaches  the  equilibrium  height.  If  the  equilibrium  plume  height 
is  higher  than  the  inversion,  the  plume  height  is  adjusted  down  to  the  inversion  height.  If  the 
stack  height  is  higher  than  the  inversion,  the  program  terminates  because  the  input 
meteorological  conditions  most  likely  do  not  apply  above  the  inversion.  If  molecular  weight  is 
not  available,  concentrations  will  be  in  mg/m3. 


3.3.12  CONCENTRATION  AVERAGING  TIME 

The  user  must  specify  the  concentration  averaging  time.  For  continuous  releases,  the 
default  value  is  15  min,  which  corresponds  to  the  short  term  exposure  limit  (STEL)  which  is 
defined  as  a  15-min.  time-weighted  average  exposure.  For  releases  of  less  than  15-nun. 
duration,  the  default  averaging  time  is  equal  to  the  release  time.  For  instantanteous  gas 
releases,  the  averaging  time  is  1  min.  The  user  can  not  enter  averaging  times  less  than  1  min. 
Also,  the  averaging  time  cannot  be  greater  than  the  release  time  since  changes  in  the 
concentrations  during  the  averaging  period  are  not  taken  into  account.  The  averaging  time 
affects  the  dispersion  coefficients  such  that  the  longer  the  averaging  time,  the  greater  the 
dispersion  coefficients,  thus  resulting  in  shorter  and  wider  plumes.  The  averaging  time  and 
dispersion  coefflcents  are  related  by  the  1/5  power  law  as  described  in  Section  2.2.4. 

3.3.13  ELAPSED  TIME  SINCE  START  OF  SPILL 

For  Instantaneous  and  finite  continuous  releases,  the  user  must  specify  the  elapsed  time 
since  th^  start  of  the  spill.  For  an  ongoing  continuous  release,  the  model  deiaults  to  a 
sufficiently  large  elapsed  time  to  assure  a  steady  state  condition  (maximum  hazard  distance). 


38 


For  an  Instantaneous  liquid  release  or  a  finite  continuous  release,  the  default  elapsed  time  is 
equal  to  the  release  time.  Except  for  short  duration  spills  (a  few  minutes),  the  default  time 
would  normally  give  the  greatest  hazard  distance.  For  instantaneous  gas  releases,  the  default 
elapsed  time  is  arbitrarily  set  at  10  min.  For  small  releases,  the  plume  may  dispersed  within 
10  min,  in  which  case  the  user  may  wish  to  enter  a  shorter  elapsed  time. 

3.3.14  TYPES  OF  OUTPUT 

The  user  can  specify  one  of  three  types  of  output. 


1 .  Toxic  corridor  plot 

2.  Concentration  at  specified  location  and  time 

3.  Maximum  concentration  at  given  height  and  time 


The  one  exception  is  when  the  no  graphics  option  is  chosen  in  SETUP.  Then  the  corridor 
plot  is  replaced  with  a  printout  of  the  hazard  distance  for  a  given  concentration(s). 

Toxic  Corridor  Plot 

If  the  user  chooses  not  to  use  the  default  concentration,  he  may  specify  up  to  three  contours 
in  any  order.  The  concentrations  can  be  specified  in  either  mg/m3  or  ppm.  The  exception  is 
that  when  the  molecular  weight  Is  not  known  the  concentration  is  in  mg/m3.  The  default 
concentration  ia  the  short  term  exposure  limit  (STEL)  in  ppm.  if  available,  otherwise  the  time 
weighted  average  (TWA)  exposure  limit.  If  neither  is  available  then  the  user  must  enter  one  or 
more  concentrations.  The  user  must  also  input  the  height  of  interest  or  choose  the  default 
height  of  2  m  (6  ft). 

AFTOX  computes  the  Y  position  of  the  concentration  of  interest  for  increasing  values  of  X. 
X  varies  in  increments  of  1 00  to  400  m  depending  on  how  rapidly  Y  Is  varying  in  the  X 
direction.  For  slowly  varying  Y.  X  Increments  approach  400  m.  The  model  starts  the  Y 
computation  at  100  m  from  the  source.  However,  if  the  centerline  concentration  at  100  m  is 
less  than  the  specified  concentration  of  Interest,  calculations  start  at  10  m  and  progress 
outward  in  10  m  increments.  Plumes  <  10  m  in  length  will  not  be  plotted.  The  model  first 
computes  the  concentration  at  the  centerline  and  then  moves  outward  in  the  Y  direction  in 
20  m  Increments  when  X>100  m,  and  2  m  Increments  when  X<100  m.  The  point  at  which  the 
computed  concentration  Is  less  than  the  concentration  of  interest  defines  the  Y  position.  The 
X  Y  positions  are  stored  in  CONCXY.DAT. 

The  model  first  computes  the  appropriate  scale  by  calculating  centerline  concentrations  for 
the  specified  time  and  height.  The  concentration  contours  arc  then  plotted  on  the  screen  as 
computations  are  taking  place.  The  contours  are  plotted  in  order  of  increasing  concentrations 
no  matter  in  what  order  they  are  entered.  These  contours  represent  average  distances. 
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meaning  that  50  percent  of  the  time  the  actual  distance  may  be  greater  than  shown  and  50 
percent  of  the  time  may  be  less  than  shown.  The  90  percent  hazard  area  for  the  lowest 
concentration  is  outlined  on  the  contour  plot  and  represents  the  area  within  which  the  plume 
is  confined  90  percent  of  the  time.  This  90  percent  area  represents  the  toxic  corridor. 

Examples  of  the  toxic  corridor  plots  are  given  In  Appendix  B. 

The  plume  plot  represents  the  position  and  size  of  a  plume  at  the  specified  tune  after  the 
start  of  the  release.  Depending  on  the  time  after  release,  this  may  not  be  the  maximum 
distance  that  the  plume  will  extend  downwind,  expedally  if  It  is  an  instantaneous  release.  The 
90  percent  hazard  area,  however,  represents  the  area  at  the  time  when  the  plume  reaches  its 
maximum  distance  downwind.  It  may  or  may  not  be  at  the  specified  time  after  release.  Run  3 
in  Appendix  B  is  an  example  of  the  maximum  distance  occurring  after  the  specified  time  after 
release.  For  more  discussion  on  the  confidence  limits,  the  reader  should  refer  to  Section 
3.3.15. 

Once  the  contour  plot  is  completed,  the  user  may  proceed  to  make  any  of  the  changes  listed 
below,  obtain  a  printer  plot,  run  another  case,  or  terminate  the  program. 


1.  time  and  meterologlcal  conditions 

2.  source  conditions 

3.  concentration  averaging  time 

4.  elapsed  time  since  start  of  spill 

5.  concentration  contours 

6.  height  of  interest 

7.  scale 


Concentration  at  a  Specified  Location  and  Time 

If  this  option  is  chosen,  the  user  must  enter  the  downwind  and  crosswind  distances,  the 
height,  and  the  time  after  start  of  the  release.  The  model  then  computes  the  concentration  for 
that  particular  point  and  time.  The  user  may  then  change  the  location  and/or  time,  choose 
another  option,  run  another  case,  or  terminate  the  program. 

Maximum  Concentration  at  a  Specified  Height  and  Time 

This  option  will  give  the  user  the  maximum  concentration  and  Its  location  for  a  specified 
height  and  time.  The  user  enters  the  height  of  Interest  and  the  time  from  start  of  release.  If  it 
Is  a  continuous  spill  that  Is  still  taking  place  and  the  specified  height  is  the  release  height,  the 
maximum  concentration  will  be  at  the  source.  In  this  case,  AFTOX  computes  the 
concentration  at  30  m  from  the  source.  If  the  specified  height  is  different  than  the  release 
height,  or  If  the  time  after  start  of  release  is  greater  than  the  duration  of  the  release,  the 
maximum  concentration  will  most  likely  occur  at  a  distance  greater  than  30  m  from  the 
source.  The  location  of  the  maximum  concentration  Is  determined  within  a  +/-5  m  accuracy 
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along  the  X  axis,  and  of  course  would  be  located  on  the  centerline  of  the  plume  at  the  specified 
height.  When  completed,  the  user  may  change  the  height  and/or  time,  choose  another  option, 
run  another  case,  or  terminate  the  program. 


3.3.15  CONFIDENCE  LIMITS 

The  model  predicts  the  mean  hazard  distance.  However,  operationally,  one  would  like  to  be 
at  least  90  percent  confident  that  the  actual  hazard  distance  will  not  exceed  the  predicted 
distance.  The  model's  concentration  contour  plot  also  shows  the  90  percent  confidence  level 
hazard  area,  or  toxic  corridor.  The  evaluation  study  described  in  Kunkel7  concluded  that  the 
predicted  hazard  distance  must  be  multiplied  by  2.1  to  be  90  percent  certa,  that  the  actual 
will  not  exceed  the  predicted  distance.  Earlier  versions  of  AFTOX  showed  the  90  percent 
hazard  area  only  for  continuous  releases.  In  AFTOX  4.0,  the  90  percent  hazard  area  Is  shown 
for  all  spills  even  though  the  2.1  factor  was  derived  from  continuous  release  data.  It  should  be 
pointed  out  that  in  deriving  the  90  percent  hazard  area,  it  is  assumed  that  the  input  data  is 
correct.  If  errors  In  the  Input  data  were  to  be  taken  Into  account,  the  toxic  corridor  would  be 
considerably  longer.  Zettlemoyer26  examined  the  effect  that  data  input  uncertainties  have  on 
the  concentration  uncertainties  using  a  Monte  Carlo  simulation  technique.  He  looked  at  wind 
speed,  emission  rate,  spill  height,  and  the  horizontal  and  vertical  dispersion  coefficients.  He 
concluded  that  errors  in  the  wind  speed  had  greater  effect  on  the  concentration  uncertainty 
than  the  other  parameters.  The  greatest  uncertainties  occurred  within  one  kilometer  of  the 
source. 

The  method  used  to  determine  the  width  of  the  hazard  area  for  the  90  percent  confidence 
level  remains  similar  to  the  method  used  by  the  Air  Weather  Service  (see  Kahler  et  al27  ).  If  the 
measured  wlrd  Is  less  than  1.8  m/sec,  (<3.5  kt),  the  hazard  area  is  a  circle  of  radius  equal  to 
2.1  times  the  predicted  hazard  distance.  If  the  standard  deviation  of  wind  direction  (o#)  Is 

known  and  the  wind  speed  is  greater  than  1 .8  m/sec,  then  the  toxic  corridor  width  (W)  is  equal 
to: 


W  =  6oe. 


(62) 


The  measured  oe  is  adjusted  by  the  one-fifth  power  law  to  the  time  duration  of  the  spill,  up  to  a 

maximum  of  one  hour.  If  the  width  is  calculated  to  be  less  than  30  degrees.  It  is  set  at  30 
degrees. 

If  o0  Is  not  known,  then  the  following  rules  apply: 


28  Zettlemoyer,  M.D.  (1990)  An  Attempy  to  Estimate  Measurement  Uncertainty  in  the  Air 
Force  Toxic  Chemical  Dispersion  (AFTOX)  Model  Master  of  Science  Thesis,  The  Florida  State 
University. 

27  Kahler,  J.P.,  Curry,  R.G.,  and  handler,  RA.  (1980)  Calculating  Toxic  Corridors. 
AWS/TR-80/003.  AD  Al 01267. 
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1)  For  neutral  or  stable  conditions,  that  Is,  a  stability  parameter  of  3.0  or  greater,  the  width 
Is  equal  to  90  degrees  for  measured  winds  of  1.8  to  9.15  m/sec  (3.5  to  10  kt).  For  winds  greater 
than  5.15  m/sec  (10  kt),  the  width  is  equal  to  45  degrees. 

2)  For  unstable  conditions,  the  width  is  equivalent  to  that  derived  from  the  ot  for  the 

different  stabilities  shown  in  Table  3.  Equation  (63)  shows  the  relationship  between  the  width 
(W)  and  the  stability  parameter  (STB). 


W«  165  -  30  STB. 


(63) 


The  width  will  vary  from  00  degrees  for  neutral  conditions  to  1 50  degrees  for  very  unstable 
conditions. 
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Appendix  A 


Chemical  Data  Base 


9 


4 


> 


1  Aerozine-50 

2  Hydrazine 

3  Hydrazine  (54%) 

4  Hydrazine  (70%) 

5  Monomethylhydrazine  (MMH) 

6  Dlmethylhydrazlne  (UDMH) 

7  JP-9 

8  JP-10 

9  Fuming  Nitric  Acfa  (IRFNA) 

10  Nitrogen  Tetroxide 

1 1  Acetaldehyde 

12  Acetone 

13  Acetonitrile 

14  Acrolein  (2-propenal) 

15  Acrylic  Acid 

16  Allyl  Alcohol  (2-propene-l-Ol) 

17  Allyl  Chloride  (3-Chloro-l-Propene) 

18  Ammonia 

19  Ammonia  (29%) 

20  Aniline 

21  Benzene 

22  Benzoyl  Chloride 

23  Bromine  Pentafluorlde 

24  Bromoforrn 

25  Bromomethane 

26  1  3-Butadlene 

27  Butane 

28  n-Butyl  Alcohol  (1 -Butanol) 

29  sec-Butyl  Alcohol  (2-Butanol) 

30  t- Butyl  Alcohol 

31  Butylene 

32  Butyraldehyde 

33  Carbon  Disulfide 

34  Carbon  Monoxide 

35  Carbon  Tetrachloride 

36  Chlorine 

37  Chlorine  Pentafluorlde 

38  Chlorine  Trlfluorlde 

39  Chlorobenzene 

40  Chloroethane 

41  Chloroform 

42  Chloromethane 

43  o-Cresol 


44  m-Cresol 

45  p-Cresol 

46  Cumene  (Isopropyl  Benzene) 

47  Cyclohexane 

48  D-D  Soil  Fumigant 

49  Dlallyl  Amine 

50  Dibutyl  Phtalate 

51  o- Dichlorobenzene 

52  m- Dichlorobenzene 

53  p- Dichlorobenzene 

54  Dlchlorodimethylsllane 

55  12  Dlchloroethylene  (els) 

56  1  2- Dlchloroethylene  (trans) 

57  1  2-Dlchloropropane 

58  Dlethanol  Amine 

59  Dimethyl  Phthalate 

60  1  4-Dicpcane 

61  Di-t-Butylethyl  Diamine 

62  Eplchlorohydrln 

63  Ethanol 

64  Ethyl  Acrylate 

65  Ethyl  Benzene 

66  Ethylene 

67  Ethyl  Chloride 

68  Ethylene  Dibromlde 

69  Ethylene  Dichloride 

70  Ethylene  Glycol 

71  Ethylene  Oxide 

72  2-Ethyl  Hexanol 

73  Fluorine 

74  Formaldehyde  (Pure) 

75  Formaldehyde  (37%) 

76  Formaldehyde  (56%) 

77  Freon  12 

78  Hydrogen  Chloride 

79  Hydrogen  Cyanide 

80  Hydrogen  Fluoride 

81  Hydrogen  Sulfide 

82  Isoamylene  (2-Methyl-2-Butene) 

83  Isobutane  (2- Methyl  propane) 

84  Isobutyraldehyde 

85  Isoprene  (2-Methyl- 1  3-Butadiene) 

86  Isopropanol  (2-Propanol) 


47 


87  Isopropyl  Ether 

88  Mesityl  Oxide 

89  Methane 

90  Methanol 

91  Methyl  Iodide 

92  Methylene  Chloride 

93  Methylethyl  Ketone 

94  Methyl  Isobutyl  Ketone 

95  Methyl  Methacrylate 

96  Nitric  Add  (Pure) 

97  Nitrobenzene 

98  Nitrogen  Dioxide 

99  Nitrogen  Trlfluorlde 

100  Nltropropane 

101  N  N-Dlmethylanlllne 

102  Oxygen  Difluorlde 

103  Perchloryl  Fluoride 

104  Pentaborane 

105  Phenol 

106  Phosgene 

107  Propane 

108  Proplonaldehyde 

109  Propylene 

110  Propylene  Oxide 

111  Pyridine 

1 12  Styrene 

113  Sulfur  Dioxide 

114  112  2-Tetrachloroethane 

1 15  Tetrachloroethylene 

116  Titanium  Tetrachloride 

117  Toluene 

1 18  o-Toluldlne 

i  19  11  1  -Triehloroethf  me 

120  1  1  2-Trichloroetha  ru- 

121  Trichloroethylene 

122  Trichlorotrifluoroethiu'e 

123  1  2  4-Trimethylbenzene 

124  Vinyl  Acetate 

125  Vinyl  Chloride 

126  Vlnylldene  Chloride 

127  o-Xylene 

128  m-Xylene 

129  p-Xylene 
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Example  Model  Runs 


USAF  Toxic  Chemical  Dispersion  Model 
AFTOX 

Hanscom  AF3  MA 
Date:  02-26-1391 
Time:  1028  LSI 

Continuous  release 
Hydrogen  Chloride 

Time  Weighted  Average  (TWA)  is  not  available 
Short  Term  Exposure  Limit  (STEL)  Is  3  ppm  (5  MG  M-3) 

Temperature  =  20  C 

Wind  direction  *  240 

Wind  speed  *  4  M/S 

Sim  elevation  angle  is  34  degrees 

Cloud  cover  Is  5  eighths 

Cloud  type  is  high  (Cl.  Cc,  Cs) 

Ground  is  wet 

T  le  Inversion  height  Is  300  M 
Atmospheric  stability  parameter  is  3.14 
Spill  site  roughness  length  is  10  cm 

This  is  a  gas  release 
Spill  Is  at  ground  level 
Emission  rate  Is  100  KG/mhi 
Elapsed  time  of  spill  Is  60  min 
Total  amount  spilled  Is  6000  KG 
Concentration  averaging  time  Is  16  min 
Elapsed  time  since  8 tan.  of  spill  is  60  min 
Height  of  interest  is  2  M 


At  60  min.  the  maximum  distance  for  5  PPM  is  1 .53  KM 
Maximum  toxic  corridor  length  =  3.24  KM  at  60  min 
Direction  &  w*ith  60  +/-  35  deg 
At  60  min,  the  maximum  distance  for  15  PPM  is  812  M 


J  CONG  MAX  DiST 

6  °PM  153  KM 
15  PPM  812M 
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USAF  Toxic  Chemical  Dispersion  Model 
AFTOX 


Hanscom  AFB  MA 
Date:  02-26-1991 
Time:  1  HOLST 
Continuous  release 
Aerozine*50 

Short  Term  Exposure  limit  (STEL)  Is  .48  PPM  (1.15  MG  M-3) 

Time  Weighted  Average  (TWA)  Is  .06  PPM  (.14  MG  M-3) 

Temperature  =  50  F 

Wind  direction  *  240 

Wind  speed  *  8  knots 

Sun  elevation  angle  is  37  degrees 

Cloud  cover  Is  2  eighths 

Cloud  type  is  high  (Cl.  Cc,  Cs) 

Ground  Is  dry 

There  Is  no  Inversion 

Atmospheric  stability  parameter  Is  2.59 

Spill  site  roughness  length  Is  50  cm 

This  Is  a  liquid  release 

Spill  is  100  gal/min 

Chemical  Is  still  leaking 

Area  of  spill  is  4000  sq  ft 

Evaporation  rate  Is  44.76  lbs/mln 

Concentration  averaging  time  Is  15  min 

Height  of  Interest  Is  6  ft 


The  maximum  distance  for  .48  PPM  is  3504  ft 
Maximum  toxic  corridor  length  =  1.39  mi 
Direction  &  width  60  +/-  43  deg 


I 
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USAF  Toxic  Chemical  Dispersion  Model 
AFTOX 


Han  scorn  AFB  MA 
Date:  02-26-1991 
Time:  1244  LST 
Instantaneous  release 
Nitrogen  Tetroxlde 

Short  Term  Exposure  Limit  (STEL)  is  2  PPM 

Time  Weighted  Average  (TWA)  is  .25  PPM 
Temperature  »  75  F 
Wind  direction  »  300 
Wind  speed  ■  7  knots 
Sun  elevation  angle  is  37  degrees 
Cloud  cover  is  0  eighths 
Ground  is  wet 
There  is  no  inversion 
Atmospheric  stability  parameter  is  3.07 
Spill  site  roughness  length  Is  10  cm 
This  is  a  gas  release 
Spill  is  at  ground  level 
Total  amount  released  is  1000  lbs 
Concentration  averaging  time  is  1  min 
Elapsed  time  since  start  of  spill  is  10  min 
Height  of  interest  is  6  ft 


The  maximum  concentration  is  7.81  PPM  (14.67  MG  M-3) 
It  occurs  at  x  *  2609  M  (8560  ft) 

Elapsed  time  since  start  of  spill  is  10  min 
Height  of  Interest  is  6  ft 


At  10  min.  the  maximum  distance  for  2  PPM  is  1.82  ml 
Maximum  toxic  corridor  length  =  6. 18  ml  at  18  min 
Direction  &  Width  120  +/-  36  deg 
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CONC  MAX  DI3T 
2  PPM  1*2  ml 


USAF  Toxic  Chemical  Dispersion  Model 
AFTOX 

Vandenberg  AFB 
Date:  02-26-1991 
Time:  1305  LST 
Instantaneous  release 
Dlmethylhydrazlne  (UDMH) 

Short  Term  Exposure  Limit  (STEL)  Is  .48  PPM  (1.15  MG  M-3) 

Time  Weighted  Average  (TWA)  is  .06  PPM  (.14  MG  M-3) 
Temperature  ■  68  F 
Wind  direction  »  270 
Wind  speed  «■  108  knots 

Standard  deviation  of  wind  direction  ■  15  deg 

Standard  deviation  averaging  time  =  15  min 

Sun  elevation  angle  is  44  degrees 

Cloud  cover  la  0  eighths 

Ground  is  dry 

There  is  no  inversion 

Atmospheric  stability  parameter  is  1.46 

Splil  site  roughness  length  is  100  cm 

This  is  a  liquid  release 

Total  amount  spilled  Is  1000  M3 

Area  of  spill  is  10000  sq  ft 

Evaporation  rate  is  56.64  KG/mln 

The  chemical  will  evaporate  in  13970.2  min 

Concentration  averaging  time  is  15  min 

Height  of  interest  Is  6  ft 


The  maximum  distance  for  .48  PPM  is  2539  ft 
Maximum  toxic  corridor  length  «  1  ml 
Direction  &  width  90  +/-  59  deg 


5  CONC  MAX  0I9T 

.48  PPM  2638  FT 


USAF  Toxic  Chemical  Dispersion  Model 
AFTOX 


Han  scorn  AFB  MA 

Date:  02-26-1991 

Time:  1313  LST 

Continuous  buoyant  plume 

Temperature  -  20  C 

Wind  direction  ■  280 

Wind  speed  =  6  M/S 

Sun  elevation  angle  is  36  degrees 

Cloud  cover  Is  0  eighths 

Ground  Is  wet 

The  Inversion  height  is  500  M 
Atmospheric  stability  parameter  Is  3.38 
Emission  rate  (KG/mln)  =  50 
Effluent  is  still  being  emitted. 

Stack  height  above  ground  (M)  *  50 

Gas  stack  temp  (C)  =  200 

Volume  flow  rate  (M3/min)  =  5000 

Effective  plume  height  (M)  =  132  at  distance  (M)  «  747 

Concentration  averaging  time  is  15  min 

Height  of  Interest  is  2  M 


The  maximum  distance  for  .2  MG/M3  Is  11.05  KM 
Maximum  toxic  corridor  length  =  23.2  KM 
Direction  &  width  100  +/-  31  deg 
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